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Abstract 

We consider the spatially homogeneous Boltzmann equation for inelastic hard 
spheres, in the framework of so-called constant normal restitution coefficients a G [0, 1] . 
In the physical regime of a small inelasticity (that is a £ [a*, 1) for some constructive 
a* > 0) we prove uniqueness of the self-similar profile for given values of the resti- 
tution coefficient a G [a*, 1), the mass and the momentum; therefore we deduce the 
uniqueness of the self-similar solution (up to a time translation) . 

Moreover, if the initial datum lies in L\, and under some smallness condition on 
(1 — a*) depending on the mass, energy and L\ norm of this initial datum, we prove 
time asymptotic convergence (with polynomial rate) of the solution towards the self- 
similar solution (the so-called homogeneous cooling state). 

These uniqueness, stability and convergence results are expressed in the self-similar 
variables and then translate into corresponding results for the original Boltzmann 
equation. The proofs are based on the identification of a suitable elastic limit rescaling, 
and the construction of a smooth path of self-similar profiles connecting to a particular 
Maxwellian equilibrium in the elastic limit, together with tools from perturbative 
theory of linear operators. Some universal quantities, such as the "quasi-elastic self- 
similar temperature" and the rate of convergence towards self-similarity at first order 
in terms of (1 — a), are obtained from our study. 

These results provide a positive answer and a mathematical proof of the Ernst- 
Brito conjecture [16] in the case of inelastic hard spheres with small inelasticity. 
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1 Introduction and main results 



1.1 The model 

We consider the spatially homogeneous Boltzmann equation for hard spheres undergoing 
inelastic collisions with a constant normal restitution coefficient a £ [0, 1) (see [17\ [U 
|23| IM]). More precisely, the gas is described by the distribution density of particles 
/ = f t = f(t,v) > with velocity v <E R N (N > 2) at time t > and it satisfies the 
evolution equation 

(1.1) |£ = Q a (fJ) m (0,+»)xR w , 

(1.2) /(0,-) = /in in R N . 



The quadratic collision operator Q a (f, f) models the interaction of particles by means 
of inelastic binary collisions (preserving mass and momentum but dissipating kinetic en- 
ergy). We define the collision operator by its action on test functions, or observables. 
Taking ip = t/j(v) to be a suitably regular test function, we introduce the following weak 
formulation of the collision operator 

(1.3) / Q a (g,f)i>dv = b\u\g*f(i>' - ip)dadvdv*, 

where we use the shorthand notations / := f(v), g* := g(v*), ip' := ip(v'), etc. Here and 
below u = v—v* denotes the relative velocity and v', i£ denotes the possible post-collisional 
velocities (which encapsule the inelasticity of the collision operator in terms of a). They 
are defined by 

, w v! , w u' 
(1-4) v = - + -, v. = ---, 

with 



a\ ( 1 + a 



w = v + u = I — - — .„+ I — - — j \u\ a. 

We also introduce the notation x = x/\x\ for any x G M , i / 0. The function b = b(u ■ a) 
in (|1.3|) is (up to a multiplicative factor) the differential collisional cross-section. We 
assume that 

(1.5) b is Lipschitz, non-decreasing and convex on (—1, 1) 
and that 

(1.6) lb m ,b M e (0,oo) s.t. Vxe[-l,l], b rn <b(x)<b M . 
Note that the "physical" cross-section for hard spheres is given by (see [T71 [13] ) 

(1.7) b(x) = b' (l-x)- E f 1 , 6' G(0,oo), 

so that it fulfills the above hypothesis (|1.5|1.6|) when N = 3. These hypothesis are needed 
in the proof of moments estimates (see |23l Proposition 3.2] and [241 Proposition 3.1]). 
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We also define the symmetrized (or polar form of the) bilinear collisional operator Q c 
by setting 



(li 



Qa(9,h)ip dv = — \ \\ b \u\ h Aip da dv dv* , 

2 J J Jrn xR n xS n-i 

k with = (i/j' + ^-V-^*)- 



In other words, Q a (g, h) = (Q a (g, h) + Q a (h, g))/2. The formula (jl.3p suggests the natural 
splitting Q a = — Q~ between gain and loss part. The loss part Q~ can be defined in 
strong form noticing that 



(Q a (g,f), i/>) = / / / b\u\g* fijdadvdv* =: (f L(g), ip), 

J J JIR^xR^xS^- 1 

where (•, •) is the usual scalar product in L 2 and L is the convolution operator 



(1.9) L(g)(v) = (b \-\*g)(v) = b / g(y*) \v - «*| dv*, with b = / b(a\)da. 

JRN Js N ~ 1 

In particular note that L and Q~ = Q~ are indeed independent of the normal restitution 
coefficient a. 

The Boltzmann equation (jl.ip is complemented with an initial datum (|1.2p which 
satisfies 

< /in E -L 1 (R Ar ), p(/ in ) := / f in dv = pe (0, oo) 

(1.10) <| , Jrn 

fin V dv = 0, £(/in) : = / /in \v | 2 cfo < OO. 

As explained in [231 02], the operator (jl.3p preserves mass and momentum, and so 
does the evolution equation: 



(1.11) 



d_ 

dt 



J, [ I 1 ,/,--(). 



while kinetic energy is dissipated 



(1.12) 



±£(f t ) = -(l-a 2 )D £ (f t ). 



The energy dissipation functional is given by 

DeU) :=b i ff* H 3 dvdv m , 

J JR N xR N 

where b\ is (up to a multiplicative factor) the angular momentum defined by 
(1.13) bi:=- f (1 - (u • a)) b(u ■ a) da. 

8 J§N-1 



In order to establish (|1.12p we have used (jl.8p and the elementary computation 

1-a 2 



A\.\2(v,v*,a) 
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(1- {u-a))\u\\ 



The study of the Cauchy theory and the cooling process of (|l.ip - (jl.2p was done in [25] . 
The equation is well-posed for instance in L\: for < /; n G L\, there is a unique global 
solution in C{R + ;L\) n L (M + ;Lg) (see Subsection 11.51 for the notation of functional 
spaces). This solution preserves mass, momentum and has a positive and decreasing 
kinetic energy. Moreover, as time goes to infinity, it satisfies: 

(1.14) £(t)^0 and f(t, •) S v=0 in M 1 (R JV )-weak *, 

where M l (K N ) denotes the space of probability measures on M. N . 



1.2 Introduction of rescaled variables 

Let us introduce some rescaled variables (which can be found in [151 [8[ [24"] for instance), 
in order to study more precisely the asymptotic behavior (|1.14p of the solution. For any 
solution / to the Boltzmann equation (jl.lh . we may associate for any r G (0, oo) the 
self-similar rescaled solution g by the relation 

g(t,v)=e- N ^f(^±,e-^v 

Using the homogeneity property (3 a (<?(A.), g(A-))(t>) = A~( N+1 ^ Q a (g, g)(Xv), it is straight- 
forward that g satisfies the evolution equation 

(1-15) ^=Qa(j9,g)-TV v -(vg). 

Any non-negative steady state < G = G(v) of (|1.15p . that is G satisfying 

(1.16) Q a (G, G) — tV v ■ (v G) = 0, 

is called a self-similar profile. It translates into a self-similar solution (or homogeneous 
cooling state) F of the original equation (II. ip by setting 

(1.17) F(t,v) = (V + Tt) N G((V + Tt)v), 

for a given constant Vq G (0,oo). Reciprocally, let us consider a self-similar solution F of 
the original equation (jl.ip . That means a solution F of (jl.ip with the specific shape 

(1.18) F(t,v) = V(t) N G{V(t)v) 

for some given non-negative distribution G = G(v) and some C , positive, increasing 
time rescaling function V(t). One can easily show (see for instance |24[ section 1.2]) that 
V(t) = rt + Vq for some constants r, Vq > and G satisfies (|1 . 16[) associated to the 
velocity rescaling parameter r. For a given self-similar profile G, associated to a velocity 
rescaling parameter r and with mass p and energy £, we may associate a new self-similar 
profile G, associated to a velocity rescaling parameter f and with mass p by setting 

G(v) = KG(Vv), V=~^-, K = V N ^. 

pr p 

The energy of G is then £ = (5) 2 £. We thus see that there exists a two real parameters 
family of self-similar profiles which can be either parametrized by (p, r) or by (p,£)- For 
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fixed mass, changing the velocity rescaling parameter r in (|1.16p corresponds to a change 
of the energy of the profile, or equivalently to an homothetic change of variable of the 
solution. Therefore it is no restriction to choose arbitrarily this constant. Also note that 
modifying Vq just corresponds to a time translation in the self-similar solution F defined 
by dLlZD. 

It was proved in [244 Theorem 1.1] that for any inelastic parameter a G (0,1), mass 
p £ (0, oo) and (thanks to the preceding discussion) any velocity rescaling parameter 
r G (0, oo), there exists at least one positive and smooth self-similar profile G with given 
mass p and vanishing momentum: 

' Q a (G,G) - tV„ • (vG) = in mN 

( L19 ) 

Gdv = p, / Gvdv = 0, 0<GeS(R N ) 



where 5(R ) denotes the Schwartz space of C°° functions decreasing at infinity faster 
than any polynomial. 

Finally, for any solution g to the Boltzmann equation in self-similar variables (|1.15p . 
we may associate a solution / to the evolution problem (11. ip . defining / by the relation 



N i ln(V + rt) 



(1.20) f(t,v) = (V + rt) N g L } (y +rt)v 



1.3 Rescaled variables and elastic limit a — > 1 

We now make the choice 

(1.21) r = Ta = p(l-a), 

and denote by G a a solution to the problem (|1.19p . At a formal level, it is immediate that 
with this choice of scaling, in the elastic limit a — ► 1, the equation (|1.19p becomes 



(1.22) 



Qi(Gi,Gi) = in R N , 

[ Gxdv = p, [ G lV dv = 0, < Gi e S(R N ) 



Moreover, multiplying the first equation of (|1.19p by \v\ 2 , integrating in the velocity vari- 
able as in (|1.12p and taking into account the additional term coming from the additional 
drift term in (|1.15p . one gets 

(1.23) 2 (1 - a)pS{G a ) - (1 - a 2 )D £ (G a ) = 0. 

Dividing the above equation by (1 — a) and passing to the limit a — > 1, one obtains 

(1.24) p£(G 1 )-D £ (G 1 ) = 0. 

It is straightforward (see Proposition 13.61 below) that the only function satisfying the 
constraints (|1.22p and (jl.24p is the Maxwellian function 

(1.25) G 1 := M 5i = M pfitSl 
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where, for any p, 6 > 0, u £ R , the function M PtUt g denotes the Maxwellian with mass p, 
momentum u and temperature 6 given by 



and where the temperature 9\ 6 (0, oo) is given by (we recall that b\ is defined in (jl,13p ) 

(1-27) h = ^ ( I M lj0 ,i(v) \v\ 3 dv 

For instance in dimension N = 3 we obtain 

s - 9n 

^ - 646f 

Moreover, in the particular case of the hard-spheres cross-section (jl.7p in dimension 3, we 
find b\ = 6' (47r)/3 and therefore 

81 



1024tt (b' ) 2 ' 



1.4 Physical and mathematical motivation 

For a detailed physical introduction to granular gases we refer to [S]. As can be seen from 
the references included in the latter, granular flows have become a subject of physical 
research on their own in the last decades, and for certain regimes of dilute and rapid flows 
this studies are based on kinetic theory. By contrast, the mathematical kinetic theory 
of granular gas is rather young and began in the late 1990 decade. We refer to [23|. [24] 
for some (short) mathematical introduction to this theory and a (non exhaustive) list 
of references. As explained in these papers, granular gases are composed of grains of 
macroscopic size with contact collisional interactions, when one does not consider other 
additional possible self-interaction mechanisms such as gravitation - for cosmic clouds for 
instance - or electromagnetism - for "dusty plasmas" for instance -. Therefore the natural 
assumption about the binary interaction between grains is that of inelastic hard spheres, 
with no loss of "tangential relative velocity" (according to the impact direction) and a 
loss in "normal relative velocity" quantified in some (normal) restitution coefficient. The 
latter is either assumed to be constant as a first approximation (as in this paper) or can be 
more intricate: for instance it is a function of the modulus \v' — v\ of the normal relative 
velocity in the case of "visco-elastic hard spheres" for instance (see [9j), which shall be 
studied in a forthcoming work [25] . 

Simplified Boltzmann models like inelastic Maxwell molecules or pseudo inelastic hard 
spheres have been proposed (see [5]) for which existence, uniqueness and global stability of 
a self-similar profile has been shown (see [3] ) , see also [2] for similar results in the driven 
case of a thermal bath. However these models do not capture some crucial physical features 
of the cooling process of granular gas, like the tail behavior of the velocity distribution of 
the rate of decay of temperature (the so-called Haff's law). For (spatially homogeneous) 
inelastic hard spheres Boltzmann models, the existing mathematical works are: 
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• the paper [8] which shows a priori polynomial and exponential moments bounds 
on any possible self-similar profile (resp. stationary solutions), whose existence is 
assumed, for freely cooling (resp. driven by a thermal bath) inelastic hard spheres 
with constant restitution coefficient; 

• the paper [T7] which shows existence of stationary solutions for inelastic hard spheres 
driven by a thermal bath, and improves the estimates on their tails of the previous 
paper into pointwise ones in this case; 

• the paper [23] which provides a Cauchy theory for freely cooling inelastic hard spheres 
with a broad family of collision kernels (including in particular restitution coefficients 
possibly depending on the relative velocity and/or the temperature), and studies the 
question of cooling in finite time or not for these various interactions; 

• the paper [24J which shows, for freely cooling inelastic hard spheres with constant 
restitution coefficient, existence of self-similar profile(s) as well as propagation of 
regularity and damping with time of singularity. 

In this paper we want to study the self-similarity properties of Boltzmann equation for 
inelastic hard spheres. Therefore as a natural first step we consider constant restitution 
coefficient a in order to have a self-similar scaling, which translates the study of self-similar 
solutions (often called homogeneous cooling states) to the study of stationary solutions for 
a rescaled equation. We also reduce to the case of restitution coefficients a close to 1, that 
is, of small inelasticity. There are several physical as well as mathematical motivations for 
such a choice: 

• the first reason is related to the physical regime of the validity of kinetic theory: as 
explained in [9J Chapter 6] for instance, the more inelasticity, the more correlations 
between grains are created during the binary collisions, and therefore the molecular 
chaos assumption, which is at the basis of the valdidity of Boltzmann's theory, 
suggests weak inelasticity to be the most effective; 

• second as emphasized in [9] again, the case of restitution coefficient a close to 1 
has been widely considered in physics or mathematical physics since it allows to use 
expansions around the elastic case, and since conversely it is an interesting question 
to understand the connection of the inelastic case (dissipative at the microscopic 
level) to the elastic case ( "hamiltonian" at the microscopic level); 

• finally this case of a small inelasticity is reasonable from the viewpoint of applica- 
tions, since it applies to interstellar dust clouds in astrophysics, or sands and dusts 
in earth-bound experiments, and more generally to visco-elastic hard spheres whose 
restitution coefficient is not constant but close to 1 on the average. 

In this framework we shall show uniqueness and attractivity of self-similar solutions 
(in a suitable sense), and thus give a complete answer to the Ernst-Brito conjecture [HI 
(stated there for the simplified inelastic Maxwell model), for inelastic hard spheres with 
a small inelasticity. Moreover we give precise results about the elastic limit and deduce 
some quantitative informations about the weakly inelastic case. 
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1.5 Notation 



Throughout the paper we shall use the notation (•) = \J\ + | • | 2 . We denote, for any 
p G [1, +00], g G K and weight function to : M. N — ► M. + , the weighted Lebesgue space L p q (uo) 
by 

LgH := {/ : R N » R measurable ; ||/|| i?H < +00} , 
with, for p < +00, 



\\f\\L^)={f RN \f(v)\ P (v) Pq u(v)dv 



1/p 



and, for p = +00, 



sup \f(v)\(v)^(v). 



We shall in particular use the exponential weight functions 

(1.28) m = m s>a (v) := e~ a ^ s for a G (0, 00), s G (0, 1), 

or a smooth version m(u) := e^^^ 2 ^ with £ G C°° is a positive function such that 
C(r) = W 2 for any r > 1, with s G (0, 1). 

In the same way, the weighted Sobolev space Wq' p {uj) (k G N) is defined by the norm 



\s\<k 



i/p 



and as usual in the case p = 2we denote H^{uj) = ' 2 (cj). The weight w shall be omitted 
when it is 1. Finally, for g G Lh, with k > 0, we introduce the following notation for the 
homogeneous moment of order 2k 



m/, : 



(ff) := / g\v\ 



2 k 



dv, 



and we also denote by p{g) = mo (5) the mass of g, £(g) = mi (g) the energy of g and 
by #(<?) = £(g)/(p(g) N) the temperature associated to g (when the distribution g has 
mean). For any p, £ G (0, 00), u G l w we then introduce the subsets of L l of functions of 
given mass, mean velocity and energy 

Cp,u '■= {h G L\; / hdv = p, / hvdv = pu}, 
JR N Jr n 

Cp,u,£ '■= {h G L\; I hdv = p, / hvdv = pu, / h \v\ 2 dv = £}. 
Jr n Jr n ' Jr n 

For any (smooth version of) exponential weight function m we introduce the Banach space 

L 1 (m~ 1 ) = L^m- 1 ) nC ,o- 
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1.6 Main results in self-similar variables 

Our main result, that we state now, deals with the evolution equation in self-similar 
variables 

(1-29) — =Q a {g,g) -r a V v -(vg), g(0, .) = g- m G C P)0 

and the associated stationary equation, namely the self-similar profile equation 

(1.30) Q a (G,G)-T a V v -(vG) =0, G £ C pfi . 

Theorem 1.1 There is some constructive a* £ (0, 1) such that for a £ [a*, 1], and any 
given mass p G (0, oo), we have: 

(i) For any r > 0, the equation (1.15\) admits a unique non-negative stationary solution 
with mass p and vanishing momentum. We denote by G a the self-similar profile 
obtained by fixing r = r Q (defined by $1.21]) ). 

(ii) Let define G\ = M p qq 1 the Maxwellian distribution with mass p, momentum and 
"quasi-elastic self-similar temperature" Q\ defined in (1.27) . The path of self- similar 



profiles a — > G a parametrized by the normal restitution coefficient is C 1 from [a*, 1] 
into W k ' 1 Pi L 1 (e a ' 1 '') for any k £ N and some a £ (0, oo). 

(in) For any a € [a*, 1], the linearized collision operator 

(1.31) h ^ C a h := 2 Q a (G a , h) — T a V v • (v h) 

is well-defined and closed on L 1 (m~ 1 ) for any exponential weight function m with 
exponent s £ (0, 1) (defined in (.1.281) ). Its spectrum decomposes between a part which 
lies in the half-plane {Re £ < p} for some constructive p < 0, and some remaining 
discrete eigenvalue p a . This eigenvalue is real negative and satisfies 

(1.32) p a = -p{\- a) + 0{l- a) 2 when a->l. 

The associated eigenspace is of dimension 1 and then denoting by (f) a = (p a (v) the 
unique associated eigenfunction such that ||^> a ||i,i = 1 an d 4> a (0) < 0, there holds 
<j) a £ S(M. N ) (with bounds of regularity independent of a) and 

(1.33) (p a -» 0i := c (M 2 - JV0i)Gi as a-»l, 



where cq is the positive constant such that ||0i||^i = 1. Finally one has constructive 
decay estimates on the semigroup associated to this spectral decomposition in this 
Banach space (see the key Theorem \5.2\ and the following point). 

(iv) The self- similar profile G a is globally attractive on bounded subsets of L\ under some 
smallness condition on the inelasticity in the following sense. For any p,£q,Mo £ 
(0, oo) there exists a** G (a*, 1), C* G (0, oo) and rj £ (0, 1), such that for any initial 
datum satisfying 

< gm g L\ n C p>0jSo , llfifoilLi < M , 
the solution g to (1.29) satisfies 

(1.34) lift -Gall i!<e N "» f . 
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(v) Moreover, under smoothness condition on the initial datum one may prove a more 
precise asymptotic decomposition, and construct Liapunov functional for the equa- 
tion il.29\) . More precisely, there exists k* E N and, for any exponential weight 
m as defined in $1.28]) and any p,£o,Mq E (0, oo), there exists a** E and 
a constructive functional H : H k * n L 1 (m _1 ) — > K such that, first, for any initial 
datum < g- m E H k * n L 1 (m~ 1 ) n Cp,o,£o satisfying 

IbinllH^nLMm- 1 ) < M)) 

i/ie solution g to hl.29\) satisfies 

(1.35) ^(i, •) = G a + c a (t) + r a (t, •), 
with c a (t) E R and r a (i, •) E L^M ) swc/t i/iai 

(1.36) |c a (i)| <C*e*>\ \\r a (t,-)\\ L i <C*e (3/2) ^*. 
And second when the initial datum satisfies additionally 

g- m >M^e- M ^\\ 

the solution satisfies also 

1 1 — ► 7i{g{t, ■)) is strictly decreasing 

(up to reach the stationary state G a ). 

Remarks 1.2 1) All the constants appearing in this theorem are contructive, which means 
that they can be made explicit, and in particular that the proof does not use any compact- 
ness argument. Unless otherwise mentioned, these constants will depend on b, on the 
dimension N, and on some bounds on the initial datum but never on the inelasticity pa- 
rameter a E (0, 1] . 

2) Theorem \1.1\ establishes that conjectures 1 and 2 in \24[ Section 5] holds true at 
least for weak inelastic model (that means for a close enough to 1). 

3) In point (iv), the condition on the restitution coefficient depends on the mass, tem- 
perature and L\ norm of the initial distribution, but this dependence is not a perturbative 
condition of closeness to the self-similar profile. This fact relies on the so-called "entropy- 
entropy production" estimates which yields "overlinear" Gronwall-type estimates, and the 
decoupling of the timescales of energy dissipation and entropy production. 

4) In tl.36\) one can prove \\r a (t, OIIl* — C( 1 f° r o,ny £ E (1, 2). Remark that here 

we do not have the decay rate e xt on the remaining part when one "removes " from gt — G a 
the projection on the energy eigenvalue, where A < would be some constant independent 
of a related to the second non-zero eigenvalue of C a . This is due to the coupling effect of 
the bilinear term, which mixes the different part of the spectral decomposition. 

5) As a subproduct the above result provides an alternative argument to the one of \2J\ 
Section 3] to show uniform (in time and inelasticity parameter) non- concentration bounds 
on the rescaled equation, in the case of a close to 1 and a general initial datum gi n E L\ 
(whereas the proof of \24\ Section 3] was valid for all a E (0, 1) but for some initial datum 
g in £LlnLP,p£ (l,oo]J. 
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6 ) Our results show that no bifurcation occurs for the self-similar profile for a close to 
1. We do not know at now if some bifurcations occur for other values of the inelasticity 
parameter. Therefore we do not know if there is a continuous branch of self-similar profiles 
parametrized by a £ [0, 1] (even if we know from [2$ that self-similar profiles exist for 
all values of the inelasticity paramaters). The best one could say in terms of "connectiv- 
ity" from the estimates we have proved on the profile together with the classical theory of 
topological degree (see J2~3\j for instance) is that there is a set K C [0,1] x T (where T is 
for instance the set of positive functions in the Schwartz space with given mass) which is 
compact, connected, and such that for any a G [0, 1], the intersection K n {a} x T is not 
empty. 

1.7 Coming back to the original equation 

When coming back to the original equation (jl.ip with the help of (11.170 and (|1.20p . 
Theorem 11.11 translates into the 

Theorem 1.3 There is a constructive a* G (0, 1) such that for a G [a*, 1], and any given 
mass p G (0, oo), we have 

(i) Up to a translation of time there exists a unique self-similar solution F a of the 
equation U.l\) with mass p, and it is given by 

F a (t,v) = (1 + r a t) N G a ((l + T a t)v), r a = p(l - a), 

where G a was obtained in Theorem More precisely, if F a is a solution of (QHJ) 
of the form i!.18\) and of mass p, there exists t$ G R such that F a (t,v) = F a (t+to,v) 
for any t > max{0, —to} an d any v G M. N . 

(ii) The self-similar solution F a is globally attractive on bounded subsets of L\ under 
some smallness condition on (1 — a*) in the following sense. For any p,£q,Mq G 
(0, oo) there exists a** G (a*,l) and rj G (0,1) such that for any q G N there is 
c q G (0, oo) such that for any initial datum satisfying 

o< f^e Lin C p>0 ,e , \\f- m \\ L x < M , 

the solution f(t,-) to satisfies 

\\f(t, •) - F a (t, .)[| L1(H9) < c q (1 + r a t)W = c q (1 + r a t )-(l-.>)-9+O(l-o0 . 

(Hi) Moreover, there exists k* G N and, for any exponential weight m as defined in \1.28\) 
and any p,£o, Mq G (0, oo), there exists a** G (a*, 1) such that, for any initial datum 
< /in G H k * n L l (m- 1 ) n C pA £ satisfying 

ll/inllifk-nLHm- 1 ) - M ° 

the solution f to hi. 29) satisfies 

(1.37) f(t, •) = F a (t, •) + c a (t) il> a (t, •) + r a (t, •) 
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where 

t/> a (t,v) = (l+T a t) N </> a ((l+T a t)v), C a (t)=C a ^(l + T-a*) 

In this expansion, the different terms have the following asymptotic behaviors (for 
any given q > 0): 

\\Fa(t, •)IIl 1 (|i>|9) = (1 + Taty q \\G a \\ L iQ v ^, 

W^ait, = (1 + T a t)~ q \\G a \\ L l(\ v \q), 

\c a (t)\ <a(i+r Q tf«/^ = a(i + T a ty 1+ °^- a \ 

3C q >0; ||f a || zl(|t , |e) < C q {l + T a t)W^»-i = C q (l + T a tr^- q+ °^- a l 

Hence the leading term in the expansion jl.37\ ) is, as expected, the self-similar solu- 
tion, and the first order correction behond self-similarity is given by the second term, 
that is the projection onto the eigenspace of the "energy eigenvalue". 

(iv) We may make more precise Haff's law on the asymptotic behavior of the granular 
temperature (see 1241) in the following way. Under the assumptions of point (in), 
the solution f = f(t,v) to satisfies 

(1.38) £(f(t.,))= ,f iGa \ 2 +o( ' ^ 



2 ' \ _i_ _ ,n3+C(1-q) 



{i + T a ty \(i + Tat ) 

(v) Under the assumptions of point (Hi) the rescaling by the square root of the energy 
familiar to physicists is rigorously justified in the sense: the solution f = f(t,v) 
to satisfies for t — ► +oo 

^(/^/^^(/t) 1 / 2 ^-^^)^^^) 1 / 2 ^ tn L 1 . 



Remark 1.4 We see from this theorem that the convergence towards the self-similar so- 
lution in indeed faster than the convergence towards the Dirac mass (hence justifying its 
interest), but also that the speed of convergence towards this self-similar solution degen- 
erates to as a — > 1 (because T a — ► when a — ► 1). This fact is surprising, since the 
self-similar solution converges towards a stationary Maxwellian distribution in the elastic 
limit, and the latter is known to be exponentially attractive for the elastic equation (see \2T$ 
for instance). As we shall see this is related to the fact that a birfurcation occurs in the 
spectrum of the linearized collision operator at a = 1 ( namely the eigenvalue corresponding 
the kinetic energy vanishes at a = 1 whereas it is non-zero for a S [a*, 1)). This remark 
may explain the fact that in the quasi- elastic limit considered - in dimension 1 - in [10], 
it is proved that the rate of relaxation towards the self-similar solution is worse than any 
polynomial. 

Proof of Theorem \1.3l Except for points (i) and (v) this theorem is an obvious translation 
of Theorem 11.11 In order to prove (i) , one first remarks that for two given self-similar 
solutions F and F, there holds 

F(t, v) = (V + A t) N G A ((V +At)v), F(t, v) = (V + A t) N GM% + At)v), 
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A\ N „ (A 



and thus from Theorem ll.il 
We deduce 



A \A 



A . \ N „ A 



with 



v4 \ i 



In order to prove (v), we introduce the function = £ (G a )V 2 /[f (/ t )V2 (i 

+ r a £)] and 

we compute 

£(/^/ 2 /(t, £(f t ) 1/2 • ) - SlGjWGaMG*) 1 * •) 
= IKt- 1 ln(l + r a t), •) - G a (£(t) 

< \\g(r- 1 ln(l + r a t), •) - G«|| Ll + - 1| ||G a [Ui + £(t) N \\G a (£(t) •) - G a [| L i. 

Using now (fL32|) . (fL35l) . (Q6j) . (OBI) and the fact that G a is bounded in W^' 1 uniformly 
in a S (a*, 1) from Theorem 1.1 (ii), we deduce 

£{f t ) N ' 2 f(t, £{f t ) 1 ' 2 ■ ) - £(G a ) N / 2 G a (£(G a )V 2 .jf^ < G (1 + r Q r 

for some constant C G (0, oo) (which depends in particular on the upper bound on 
l|G a ||jyi,i), from which (v) follows. □ 

Remark 1.5 Let us emphasize that the temperature 8± of the limit Maxwellian G\ is 
"universal" in the sense that it depends only on the collisional cross-section b (through its 
angular momentum), and not for instance on the density distribution. 

The temperature of the self-similar solution F a = F a (t,v) associated to a self-similar 
profile G a decreases like 

B{Pa{t '' )) = (l + fa-a)tr 
Hence when a is close to 1 (small inelasticity) we obtain 

0(F a (t,-)) 



(l + p(l-a)t) 2 ' 

Therefore, as soon as the self-similar solutions correctly describe the asymptotic (at least 
in the framework of point (ii) of Theorem \1.3\) . which is conjectured by physicists, generic 
solutions satisfy 
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for an inelasticity coefficient a close to 1. 

Hence we shall denote the universal quantity 0\ as a "quasi- elastic self-similar temper- 
ature". Remark that its definition as the temperature of G\ seems to depend on the choice 
of the scaling. However changing this scaling by some asymptotically equivalent one, as 
a — ► 1, would only adds a factor which would then disappear when coming back to the 
solution to the original equation U.l\) . Therefore a more "canonical" way to define this 
quasi-elastic self-similar temperature could be 



where f a denotes a generic solution with mass p to equation hl.l}) . 
1.8 Method of proof and plan of the paper 

The first main idea of our method is to consider the rescaled equations (|1.15j) and (|1.16l) 
with an inelasticity dependent anti-drift coefficient r Q which exactly "compensates" the 
loss of elasticity of the collision operator (in the sense that it compensates its loss of ki- 
netic energy). This scaling allows by some technical estimates to prove uniform bounds 
according to a for the family of self-similar profiles G a to the equation (ll.30p . The sec- 
ond main idea consists in decoupling the variations along the "energy direction" and its 
"orthogonal direction" . This decoupling makes possible to identify the limit of different 
objects as a — > 1 (among them the limit of G a ). The third main idea is to use systemati- 
cally the knowledges on the elastic limit problem, once it has been identified thanks to the 
previous arguments. In particular we use the spectral study of the linearized problem and 
the dissipation entropy-entropy inequality for the elastic problem. This allows to argue by 
perturbative method. Let us emphasize that this perturbation is singular in the classical 
sense because of the addition of a (limit vanishing) first-order derivative operator, but 
also because of the gain of one more conservative quantity at the limit (which implies in 
particular at the linearized level that the "energy eigenvalue" \x a is negative for a^l but 
converges to [i\ = in the limit a — > 0). 

In Section 2, we use the regularity properties of the collision operator in order to 
establish on the one hand that the family (G a ) is bounded in H°° n L 1 (m _1 ) uniformly 
according to the inelastic parameter a (the key argument being the use of the entropy 
functional which provides uniform lower bound on the energy of G a ) and on the other 
hand that the difference of two self-similar profiles in any strong norm may be bounded 
by the difference of these ones in weak norm (the key idea is a bootstrap argument). This 
last point shall allow to deal with the loss of derivatives and weights in the operator norms 
used in the sequel of the paper. 

In Section 3, we prove that a t— > is Holder continuous in the norm of its graph and is 
Holder differentiable in a weaker norm. As a consequence we deduce that G a — ► G\ when 
a — * 1 with explicit "Holder" rate, which (partially) proves point (ii) Theorem 11.11 The 
cornerstone of the proof is the decoupling of the variation G a — G\ between the "energy 
direction" and its "orthogonal direction". 

In Section 4, we prove uniqueness of the profile G a for small inelasticity (point (i) of 
Theorem II. ip by a variation around the implicit function theorem. We also deduce that 
a i— ► G a is differentiable at a = 1. 
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Section 5 is devoted to the study of the linearized operator C a , and we partially inspire 
from the method of [27J . We prove point (iii) of Theorem 11.11 and we end the proof of 
point (ii) of Theorem II, 11 We obtain information on the localization of the spectrum and 
we establish some decay estimates on the associated semigroup. Let us emphasize that for 
technical reasons we state our results in an L 1 framework (because mainly we are not able 
to generalize Lemma 15.81 to an I? framework), which makes the spectral analysis more 
intricate. The proof proceeds as follows (the cornerstone idea is again the decoupling of the 
variations in the "energy direction" and its "orthogonal direction"). First, we localize the 
essential spectrum in the half plan A^ = {z £ C, Jftez < /x < 0} with the help of Weyl's 
theorem, the compactness properties of C a and the "rough" (Holder type) convergence of 
Qa(G a , •) to Qf(Gi, •) in the "good" norm of the graph. Second, we localize the discrete 
spectrum lying in A^ = {z G C, $tez > p,} in the disc {z G C, \z\ < C (1 — a)}, thanks 
to estimates on the resolvent of C a . Third we establish that the spectrum T,(C a ) of C a 
satisfies S(£ a )nA^ = {fi a }, where fi a has multiplicity 1 (the proof mainly takes advantage 
of the "precise" convergence of Q^(G a ,-) to Q± (Gi,-) in "bad" norm, together with a 
regularity estimate holding on the discrete eigenspace). Last we establish the expansion 
(jl.32p using the energy equation associated to the eigenvalue [i a . The decay properties of 
the linear semigroup are then deduced from resolvent estimates and the above localization 
of the spectrum. 

Section 6 is devoted to the proof of points (iv) and (v) in Theorem 1 1 . 1 1 which is split in 
several steps. First we establish a "linearized asymptotic stability result" by decoupling the 
evolution equation (jl.29p along the "energy direction" and its "orthogonal direction" , and 
using the semigroup decay estimates and the quadratic structure of the collision operator. 
Second we establish a "non-linear stability result" by decoupling the evolution equation 
(|1.29j) , using the energy dissipation equation along the "energy direction" and the entropy 
production method on its "orthogonal direction" (let us mention that this method follows 
closely the physical idea that for small inelasticity the "molecular" timescale of thermal- 
ization of velocity distribution decouples from the "cooling" timescale of dissipation of 
energy). Third we prove the asymptotic decomposition and we exhibit a Liapunov func- 
tional for smooth initial data (point (v)) by gathering (and slightly modifying) the two 
preceding steps. Fourth and last, we prove point (iv) for general initial data, gathering 
the previous arguments with the decomposition of solutions between a smooth part and a 
small remaining part as introduced in |28j . 

2 A posteriori estimates on the self-similar profiles 

In this section we prove various a posteriori regularity and decay estimates on the self- 
similar profiles (or the differences of self-similar profiles), uniform as a — > 1, which shall 
be useful in the sequel. 

2.1 Uniform estimates on the self-similar profiles 

For any a G (0, 1) we consider Q a the set of all the self-similar profiles of the inelastic 
Boltzmann equation (jl.ip with inelasticity coefficient a, with given mass p G (0, +oo) and 
finite energy. More precisely, we define Q a as the following set of functions 

Q Q ■= \o < G G L\ satisfying ([Onjl j. 
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For some fixed oq G (0, 1), we also define 



Q — U Q6[ao,l)^a- 

The fact that for any a G (0,1), Q a is not empty was proved in [23], where a so- 
lution of (|1.30|) was built within the class of radially symmetric functions belonging to 
the Schwartz space. Here we show that any self-similar profile G a G Q belongs to the 
Schwartz space and that decay estimates, pointwise lower bound and regularity estimates 
can be made uniform according to the inelasticity coefficient a G [cto, 1). Let us emphasize 
once again that the choice of the velocity rescaling parameter r a = p (1 — a) in fll.3Qf) is 
fundamental in order to get that uniformity in the limit a — > 1. Let us also mention that 
our choice of scaling for the equation (11.301) is mass invariant, that is G with density p(G) 
satisfies the equation if and only if G/p(G) satisfies the equation with p = 1. Therefore 
all the estimates on the profiles are homogeneous in terms of the density p. 

Proposition 2.1 Let us fix ao G (0,1). There exists 01,02,03,04 G (0, 00) and, for any 
k G N, there exists G (0, 00) such that 



(2.1) VaG[« ,l),VG a Gfc, 



||Ca[ln( e »l M) — a 2> 1 1 G a \ \ fjk qr]V ) < Cfc, 

G a > a 3 e- a4 H 8 . 



We first recall the following geometrical lemma extracted (in a slightly specified form) 
from \23\ Lemma 2.3 & Lemma 4.4], that we shall use several times in the sequel. 

Lemma 2.2 For any a G (0, 1] and a G S^ -1 we define 

<j>* a = 0* )UjCT : R N - R N , v* 1 * v' 
4><x = (t>a,v*,o : M. N —* R N , v 1 — ^ t/ 

and i/ie Jacobian functions J* = det (D <fi* av a ), J a = det (-D <t>oc ;V *,a), as we// as i/ie cone 

ft 5 = 5 ,<t = {n G M N , u ■ a > 5 - 1}, 

/or any 5 G (0, 2) and a G S^" 1 . 

For any 8 G (0, 2), 0* defines a C°° -diffeomorphism from v + onto v + wni/i 
w*((5) = 1 + y/S/2 and (f) a defines a C°° -diffeomorphism from + onto + ^ Wa (<j) 
with 

rx\ 11 5 - 1 + r a 

w a (d) = 1 + 



(l + 2(5 - l)r a + rl 



1/2 



and r a = (1 + a)/(3 — a). 

Moreover, there exist C G (0, 00) such that with C$ = C/5 

(2.2) CJ 1 \v — v*\ < \4> a (v) — v*\ < 2 \v — v*\, 

(2.3) |CV) - C,V)| < C 5 \a' - a\ \v' - n 



* 1 • 



(2-4) \Ja\<C s , iJ^l^Cs, I J~ l - J-/| < Cf |a' - a| 
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on v* + Qs, uniformly with respect to the parameters a, a' G [0, 1], a G S^ -1 and v* G l w . 
Tae same estimate holds for cff a on v + f2<5. Finally, for any a, a' G [0,1], <7 G S^ -1 , 
G 1^ and t G [0, 1], iaere holds 

(2-5) ttf + il-t) & = 

for some at belonging to the segment with extremal points a and a' . The same result holds 
for 0*. 

We will also need the following elementary result in order to estimate the convolution 
operator L defined in (jl.9p . 



Lemma 2.3 For any function g G L^M^) there exists some constants c\,c 2 G (0, oo) 
such that 

(2.6) c 1 (l + \v\)<L(g)<c 2 {l + \v\). 

Moreover, if g satisfies £{g) >a\p and m 3/ / 2 (g) < a 2 p, for some constants ai,a 2 > 0, we 
can take c\ = C _1 p, c 2 = C p in \2. 6\) for some explicit constant C > depending only 
on ai,a 2 > 0. 

Proof of Lemma [23 The upper bound in (|2.6p is immediate. As for the lower bound, we 
have, on the one hand, by Jensen's inequality, 



(2.7) / g* \u\ dv* > p\v\. 

Jr n 

On the other hand, by triangular inequality, 

/ 5* \u\ dt>* > mi/2 — \v\ m . 

Jrn 

By Holder's inequality we have > £ 2 m^ > Cq p for some explicit constant Co > 
depending only on oj, a 2 . As a consequence 

(2.8) / g*\u\dv*>p(C -\v\). 

These two lower bounds (|2.7l 12. 8|) imply immediately that 

/ g* \u\ dv* > C^ 1 p(l + \v\). 
Jrn 

for some explicit constant C > depending only on Co- □ 

Proof of Proposition \2.1l We split the proof into several steps. In Steps 1, 2 and 3, we 
establish the smoothness for any profile G a G Q as well as upper and lower bounds on its 
tail. In Steps 4, 5, 6, 7, 8 and 9, we show that these estimates actually are uniform with 
respect to the choice of the profile G a G Q a and a G [ao, 1). Thanks to Steps 1, 2 and 3 
the computations then performed are rigorously justified. 

We fix a G [ao, 1) and G a a solution of (|1.3(J|) for which we will establish the announced 
bounds. From now we omit the subscript "a" when no confusion is possible. 
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Step 1. Moment bounds. From \14\ Proposition 3.1], by taking = G in the evolution 
equation (jl . 15H . we get that G 6 L\ for any k € N. 

S'tep 2. L 2 a posteriori bound. We aim to prove that G 6 I 2 . Let us fix A > and let 
us introduce the C 1 function 

x 2 ( A 2 \ 

Aa(sc) := — 1 x <a + f A x - — J 1 x> a. 

We multiply the equation (|1.30p by h! A {G) = min{G, ^4} := Ta{G). Once again we 
shall omit the subscript "A" when no confusion is possible. After some straightforward 
computation we get 

/ (t(G)GL(G) + p(l- a)NT{G) 2 /2) dv = [ T(G)Q + (G,G)dv. 

Since L(G) > c\ (1 + \v\) thanks to Lemma and k(G) < GT(G) we have 

(2.9) ci / A(G)(1 + M)d« < / T(G)GL(G)dv 
Jr n Jr n 



< [ T(G)Q + (G,G)dv <h + h + h + h, 
Jr n 



where the terms Ik are defined in the following way, splitting the collision kernel into 
some smooth and non-smooth parts. Let O : R — ► R+ be an even C°° function such that 
support C (—1, 1), and L0 = 1. Let G : R^ — > R+ be a radial C°° function such that 
support C -B(0, 1) and J RN = 1. Introduce the regularizing sequences 



@ m (z) = mO(mz), zeR, @ n {x) = n N G(nx), xeR N . 

As a convention, we shall use subscripts S for "smooth" and R for "remainder". We 
denote $(u) := First, we set 

$S,n = 0n*($ 1a), $fl,n = * - *S,n, 

where A n stands for the annulus A n = {x £ R^ ; ^ < |x| < n}. Similarly, we set 

&S,m(2) = 6m * (& lx m ) 0), b R>m = 6 - 6s,, 



'.mi 



where I m stands for the interval Z m = {i£R; — 1 + ^ < |a?| < 1 — ^} (6 is understood 
as a function defined on R with compact support in [—1,1]). We then define 

h= I T(G)Q+(G,G)dv, 
where is the gain term associated to the cross-section Br :=\u\ &R, m , 

h= I T(G)Q+ s (G,G)dv, 
where Q~r S is the gain term associated to the cross-section Brs '■= &R, n bs, m , 



T{G) Qt(x(G), G) + Q + S {T{G), X (G) 



dv, 
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where Qg is the gain term associated to the smooth cross-section B$ := &s,n,bs,m and 
X (G) := G - T(G) and finally 



h = I T(G)Q+(T(G),T(G))dv. 
it" 



We estimate each term separately. We omit the subscripts m and n when there is no 
confusion. For I\ we proceed along the line of the proof of the estimate for the term I r in 
Proof of Theorem 2.1]. Using Young's inequality xT{y) < A(x) + A(y) we have 



h = / / / GG*T(G')b^ m \u\dv dv*da 

< /// G[A{G*) + A(G')]b Rtm l & . t7 < \u\dvdv*da 

+ / / / G* [A(G) + A(G')] b R>m l 6 . CT >o |u| cfo cfe* da = J^i + ... + J 1>4 . 

J ■/ iR^xR^xS*- 1 

We just deal with the term 1^2, the others may be handled in a similar (or even simpler) 
way. Making the change of variables — ► v' = <^j(f*) (for some fixed v, a) and using the 
elementary inequality \u\ < 4 \ v' — v\ valid when a ■ u < 0, there holds 

h,2 = \\\ G A(G') b^m l«-o-<o \ u \ dv dv* da 



= 2 N+2 GA(G')b R , m l & . a < \v-v'\dv'dvda 

J J iR^xR^xS*- 1 

< 2 Ar+2 ||^, m || Ll ||G|| L i f A(G)(l + \v\)dv. 

1 Jr n 

Since the same estimates hold for all the terms In-, we obtain 

(2.10) h <e{m)\\G\\ L i / A(G) (v) dv with e(m) — ► 0. 

i J RN m->oc 

For I2 we proceed along the line of the proof of the estimate for the term I in 
Proof of Proposition 2.5]. Using again Young's inequality xT(y) < A(x) + A(y) and the 
trivial estimate &R, n < C rT x (|w| 2 + |f*| 2 ) we get 

h = fit GG*T{G')b s , m $R,ndvdv*da 

< - III G\v\ 2 [A(G*) +A{G')}b Sm dvdv*da 

n J J JR N xR N xS N - 1 

+- / / / G*\v*\ 2 [A(G)+A(G')]bs,mdvdv*da = I 2 ,i + ... + l2A- 

n J J Jrn xR n xS n-i 

Because of the truncation on b of frontal and grazing collisions, both changes of variables 
v — ► v' = 4> a {v) (for fixed cr) and — > v' = (f° r fixed u, <r) are allowed (and the 

jacobian of their inverse is bounded). Hence in a similar way as for the term I\ we obtain 

(2.11) i 2 <^M\\ G \\ f A(G) dv. 

n 2 Jrn 
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For ^3, using again Young's inequality, plus T(G) < G and the fact that both changes 
of variables v — > v' = <p a {v) (for fixed u*,cr) and v* — ► v' = 4>* a {v*) (for fixed v,a) are 
allowed, we have 



h = /// (r(G , )+T(G'J)[Gx(G*) + x(G ! )r(G f *)]6s,m^s,nrfwdw*d<T 

J J JR N xR N xS N - 1 

< C{n) If! {x(G,)[A(G') + A(G'J + 2A(G)] 

+ X (G) [A(G') + A(Gt) + 2 A ((?*)] } 6 5 , m <fo cfc* tfV. 

We deduce as before 

(2.12) h<C m>n \\ X (G)\\ L i [ A(G) dv 

Jr n 

for some constant C m ,n > 0. 

Finally for I4, we argue as in the proof of |24| Proposition 2.6] for the treatment of the 
term involving Qg, and we get for some 9 £ (0, 1) 

(2-13) h<C m , n \\T(G)\\l+ 2e \\T(G)\\lt 6 \ 

for some constant C m ^ n > 0. 

Gathering (f!T9|) . ([2T0j) . ([2"TTj) . ([2T2|) . (|2TT5|) and taking m, next n and finally A > 
A{G) large enough we may control the terms I\, I2 and I3 by the half of the left hand side 
term of (|2.9|) (for I3 we use that Hx^C^IIl 1 ~~ > when A — > 00). Note that the condition 
A > depends on the distribution G (by the mean of some non-concentration bound), 

but shall play no role since we shall take the limit A — * +00 in the end. We obtain 



VA>A(G), % [ AA(G)(l + \v\)dv<C b . p£{G) \\T A (G)\\ 
2 Jrn ' 



2(1-0) 
L 2 



for some constant C bp£ ^ > depending on the cross-section b and on the profile G via 
its energy. Using that Ta(G) 2 /2 < Aa(G) we deduce 



VA>A(G), ^\\T A (G)\\H<C b , pAG) 



and we then conclude that G S L 2 passing to the limit ^4 — > 00 in the preceding estimate, 
with the bound 



(2-14) \\G\\ L 2 < 



1 



ci 



Remark 2.4 Note that the 1? bound {2.1J$ only defends on the distribution G by the 
mean of the energy £(G) and the constant c\. Therefore, thanks to Lemma \2.3[ this bound 
only depends on a lower bound on the energy £{g) and an upper bound on the third moment 
1*13/2(5). 

Step 3. Smoothness and positivity. Thanks to |24^ Theorem 1.3] and [SJ Theorem 1], 
taking g m = G as an initial condition in (|1.15p we have that G belongs to the Schwartz 
space of C°° functions decreasing faster than any polynomials, and that G > a\ e~ a2 for 
some constant 01,02 > 0. 
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So far the estimates in Step 3 may be not uniform on the elasticity coefficient a £ [ato, 1) 
and on the profile G a - The aim of the following steps is to prove that they actually are 
uniform. Note however that estimates of the previous steps shall ensure that the following 
computations are rigorously justified. 



Step 4- Upper bound on the energy using the energy dissipation term. We prove that 

(2.15) Va€(0,l] £<j$P- 
From equation (|1.23p on the energy of the profile G there holds 

(2.16) (l + a)h / GG*\u\ 3 dvdv* =2p G\v\ 2 dv. 

Jr n Jr n Jr n 

From Jensen's inequality 



/ 

JR> 



and Holder's inequality 



we get 



\u\ 3 G* cfo* > p \v\ 3 , 

\ 3/2 



\v\ 3 Gdv > p~ l/2 ( [ \v\ 2 Gdv) 
\Jr n J 



(l + a)6i p 1/2 £ 3/2 <2p£ 
from which the bound (|2. 15[) follows. 

Step 5. Lower bound on the energy using the entropy. We prove 

Na 4 

(2.17) Va€(0,l] £>^-p. 

o 

Remark 2.5 The choice of scaling we have made for the evolution equation in self-similar 
variable becomes clear from this computation: it is chosen such that the energy of the self- 
similar profile does not blow up nor vanishes for a — > 1. The restriction a £ [ao, 1), 
«o > ; is then made in order to get a uniform estimate from below on the energy. 

By integrating the equation satisfied by G against log G we find 

/ Q(G, G) log Gdv - p(l - a) [ logG V« • (vG) dv = 0. 
Jr n Jr n 

Then we write the first term as in \17\ Section 1.4] to find 

\\\\ GG.flog^-^ + l) Bdvdv.da 

^ J J Jr 2N xS n - 1 V GG* GG* / 

+1 III (G'Gi-GG*) Bdvdv*da + p{l-a) f v-V v Gdv = 0. 

2 J J Jr2n xS n-i J K jv 

If we denote 

(2.18) D H , a {g) = \fff ggJ^A-iogtA -A Bdvdv*do->0, 

2 J J 7R2JV xsn-i \gg* gg* J 
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(recall that in this formula the post-collisional velocities v',v* are computed according to 
the inelastic formula (jl.4p with normal restitution coefficient a G (0, 1]), we can write 

(2.19) -D Ha (G)+ (4r-l J h ff GG m \u\dvdv*-(l-a)N p 2 = 0, 

V« / J Jr 2N 

with 62 := II^IIl 1 ! aim thus we get 

GG* \u\ dv dv* = ( N p 2 + -J—D H , a (G)) > ^-/A 

1 + a \ 1 — a J I 

On the other hand, from Cauchy-Schwarz's inequality 
G G* \u\ dv dv* < 



<( f f GG^dvdv^] '(II GGM 2 dvdv^] ' = V2p z / 2 £ 1 ' 2 , 

\J J \J J R 2N J 

and then the bound (|2,17p follows gathering the two preceding estimates. 

Step 6. Upper bound on (exponential) moments using Povzner inequality. There exists 
A, C > such that 

V«e[0,l), / G(v)e Alvl dv < C p. 

We refer to [S] where that bound is obtained as an immediate consequence of the 
following sharp moment estimates: there exists X > such that 

(2.20) Va€[0,l), m fc = / G\v\ k dv <T(k + 1/2) X k/2 p. 

It is worth noticing that in [8] the Povzner inequality used in order to get (12.2QH is uniform 
in the normal restitution coefficient a £ [0, 1] and that the factor p comes from our choice 
of the scaling variables (in which p is involved). 

Step 7. Uniform upper bound on the I? norm. From (|2.17p . (|2.2U|) and Remark 12. A\ the 
I? bound (|2.14p is uniform on a G [cto, 1) and G G Q a . 

Step 8. Smoothness. It is enough to show some uniform bounds from above and below 
on the energy together with uniform non-concentration bounds on the self-similar profiles 
in Q, in the form of upper bounds on the L 2 bounds for instance. Indeed the proofs 
of [24\ Proposition 3.1, Proposition 3.2, Proposition 3.4, Theorem 3.5 and Theorem 3.6] 
then apply straightforwardly (in these proofs we did not use the part associated with the 
anti-drift in the semigroup). Therefore the uniform bounds on the H k norms for all k > 
follows from these results. 

Step 9. Pointwise Lower bound. It is a consequence of the following lemma. □ 

Lemma 2.6 Let g G C([0, 00); L3) be a solution of the rescaled equation $1.29\) with in- 
elasticity parameter a G (0, 1) and assume that for some p > 1 and C,T G (0, 00) 

sup \\g\\ L p nL i < c. 

[0,T] J 
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(i) For any t± G (0, T) there exists a\ E (0, oo) (depending on C, p and t\ but not on T) 
such that 

(2.21) Vie [ii,T], VveR N , g{t,v) > a^ 1 e^ 1 ^ . 

(ii) If furthermore, g- in satisfies 

9in(v)>a^e- a ^ 8 , 

then \2.21\) holds with to = and some constant a\ £ (0, oo) (depending on C, p,ao but 
not on T). 

Proof of Lemm& \2. 61 We only prove (i), the proof of (ii) being similar. Let us fix t\ G (0, 1). 
We closely follow the proof of the Maxwellian lower bound for the solutions of the elastic 
Boltzmann equation (see [H] [30] ) taking advantage of some technical results established 
in its extension to the solutions of the inelastic Boltzmann equation (see [241 Theorem 
4.9]). The starting point is again the evolution equation satisfied by g written in the form 

d t g + T a vV v g + (r a N + C + C\v\)g = Q+(g, g) + (C + C \v\ - L(g)) g, 

where the last term in the right hand side term is non-negative for some well-chosen 
numerical constant C £ (0, oo) thanks to Lemma 12.31 (|2.20|> and (|2.17j) . Let us introduce 
the semigroup Ut associated to the operator r a v- V v + X(v), where X(v) := r Q N+C+C\v\, 
which action is given by 

(U t h)(v) = h{ve~ Tat ) exp f- jf \{ve~ s )ds^j . 

Thanks to the Duhamel formula, we have 

(2.22) Vt>0,Vr>0, g(t + r, .) > / U t - s Q + (g(s + r, .), g(s + r, .)) ds. 

Jo 

Noticing that 

(- J \(ve- s )ds^ > -(C\v\+r a Nt + Ct), 

and repeating the arguments of Steps 2 and 3 in the proof of [231 Theorem 4.9], we get 
that 

(2.23) Vt>r, g(t,.)> V l Bm (v) 

with t = T\ = t\/2 and some constant r) = rjx > 0, S = 5\ > 1. Let us emphasize that 
here we make use of Lemma 4.6, Lemma 4.7 and Lemma 4.8 in [24J where the constants 
exhibited in these ones are uniform in a £ [ao, 1) thanks to the uniform LP n L\ estimates 
assumed on g. 

Now, on the one hand, from [24], Lemma 4.8], there exists k S (0, oo) such that 

Qt^BiO^ABiO,!)) > Kl B(0,./5/2) 

which in turns implies 

(2.24) W8 > 0, Q+(l B( o l(5 ),l fl (o,5)) > kS-"- 1 1 B (o,VE/2S)- 
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On the other hand, there exists «' G (0, oo) such that 

(2.25) V$>0, Vs€[0,l], U s (l Bm )>Ke~ cs l Bm - 

From flgjgp with r] = rji, S = S 1} and making use of (I2T22]) . (EHHD . (|235|) . we get that 
(|2T23j) holds with 

T = T2 = n + ^, 5 = 5 2 = ^-5i and n = r] 2 = (r 2 - n) k" rjj e~ c> 5l , 

where and C depends on C and N. Iterating the argument we get that (I2.23|) 

holds with r = r fc = r fc _! + ti 2~ k = (1 - 2" fc ) h, 5 = 5 k+1 = (V5/2) h+1 and 

m+l = ( ^ tl )l+2+...+2^i^ e -C'(5 fe +25 fc _ 1+ ...+ 2'=-M 1 ) 2 -[fe+2(fc-l)+...+2^il] > 
with A := V K "*i v^?i e_C " 5l /2- In other words, using that (^"^^ > 2, we have proved 

Vi > ti, V& € N, > l B(0)2 ^5i)( u )> 

from which we easily conclude. □ 

2.2 Estimates on the difference of two self-similar profiles 

In this subsection we take advantage of the mixing effects of the collision operator in order 
to show that the L 1 norm of their difference of two self-similar profiles (corresponding to the 
same inelasticity coefficient) indeed controls the H k r\L l (m~ 1 ) norm of their difference for 
any k £ N and for some exponential weight function m, uniformly in terms of a £ [ao, 1)- 

Proposition 2.7 For any k > 0, there is m = exp(— a a G (0,oo) and Ck > such 
that for any a £ [ao, 1) and any G a ,H a G Q a there holds 

(2.26) \\H a — G a ||#-k n £i( m -i) < Cfc \\H a — Galk 1 - 

Proof of Proposition \2.7\ We proceed in three steps. It is worth mentioning that all the 
constants in the proof are uniform in terms of the normal restitution coefficient a G [ao, 1), 
as they only depend on the uniform bounds of Proposition 12. II and some uniform bounds 
on the collision kernel. 

Step 1. Control of the L 1 moments. We prove first that there exists A,C £ (0, oo) such 
that 

VaG[a ,l), f \H a -G a \e Alvl dv <C [ \H a -G a \dv. 
Jr n Jr n 

Let us consider some normal restitution coefficient a G [ao, 1) and two self-similar profiles 
G, H G Q a (here again, we omit the subscript a when there is no confusion). We denote 
D = G - H, S = G + H and ip = \v\ 2p sgn(D), p G ~ N, p > 3/2, where sgn(£>) denotes 
the sign of D. The equation for D reads 

= Q a (G,G)-Q a (H,H)-p{l-a)V v -(vD) 

(2.27) = 2Q a (D,S)-p(l-a)V v -(vD). 
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Multiplying equation (|2.27[) by tp, we get 



/ 

Jr> 



BDS* 

}IV v «jV-l 



ip* + tp' — — tp dv dv* da 
-p(l-a) / V v (vD) \v\ 2p sgn(D) dv 



< [ \u\\D\S*K p dvdv* + 2 [ \u\\D\S*\v*\ 2p dvdv* 

JR N xR N JR n xR n 

+ P [ \D\v V(\v\ 2p )dv 
Jr n 

with 

K p (v,v*):= [ {\v'\ 2p + \v'f p -\v\ 2p -\v*\ 2p )b{o -u)da. 



From [8j Corollary 3, Lemma 2], there holds 

K p (y,v*) < 7p S P - (1 - 7 P ) (\v\ 2p + M 2p ) 
where (7^)^=3/2,2,... is a decreasing sequence of real numbers such that 

(2.28) 0<7 P <minjl,-^j- 



and Tip is defined by 



k=l 



E p := (I) (\v\ 2k \v*\ 2p ' 2k + \v\ 2p ' 2k \v ' 2fc 



with k p := [(p + l)/2] is the integer part of (p + l)/2 and ( ^ J stands for the binomial 
coefficient. As a consequence, 

(1 — 73/2) / M 2p \u\ 5* |£>| (iu cfo* < 7 P / |u| |D| S* £ p cfo* 

Jr^xM^ JR n xR n 

+2 / |u| |D|5* \v*\ 2p dvdv* + 2pp [ \D\ \v\ 2p dv. 
Jr n xr n Jr n 

Using Lemma [2.31 in order to estimate L{S) from below, the inequality \u\ < \v\ + \v*\ and 
introducing the notations 

d fc := f \D\\v\ 2k dv, s k := f S\v\ 2k dv, 
Jr n Jr n 

we get, for some numerical constant C G (0,oo), 

(2-29) ^d p+1/2 < 7p5 p + (d s p+1/2 + d 1/2 s p ) +2ppd p , 

with 

<Sp : = ^ ( . J (dfc+1/2 s p _ fc + d fc Sp_fc +1 / 2 + dp_fc +1 / 2 Sfc + d p _fc Sfc +1 / 2 ) • 
fe=i ^ ' 
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From Proposition 12. 11 or more precisely (|2.20p . we know that s& < pT(k+ 1/2) x k for any 
k > 1 and for some x 6 (l,oo). By Holder's inequality, we also have 

di +i <d p+ ,df. 

Repeating the proof of [U Lemma 4], for any a > 1, there exists ^4 > such that 

S'p < A p (d + di /2 )r(op + a/2 + 1) Z p 

with 

' = j-J? 9 ^ ^ fc +!/ 2 a p-ki$k & p -k+l/2, dp-k+l/2 &ki ^p-k ^k+l/2}, 

and 

, _ dfc _ s fc 

fc : ~ (do + d 1/2 )r(afc + l/2)' ak : ~ pT(ak + 1/2)' 

We may then rewrite (12.29H as 

T{a P + 1/2)1/* ^+XA, < A7p r( g^° / 1 2 / + 1) Z p + ( ffp+1/2 + a p ) + 2pp 8 p . 
On the one hand, from (|2.28j) . there exists A' such that 

On the other hand, thanks to Stirling's formula n! ~ n n e~™ \j2-nn when n — > oo and the 
estimate (12.28p . there exists A" > such that 

(1 - 7p ) T(ap + 1/2) ^ > ^> a / 2 Vp = 3/2, 2, . . . 

Therefore, 

p a/2 § l+l/2p < p a/2-V2 Zp + + ^ ^ + 2 pp^. 

We finally obtain 

d k < x k T(ak + 1/2) (d + d 1/2 ), 

and we easily conclude as in [8, Proof of Theorem 1] or in [23} Proof of Proposition 3.2, 
Step 2). 

Step 2. Control of the L 2 norms. For k = 0, the propagation of the L 2 norm is immediate 
using the result [24, Corollary 2.3]. Indeed one just has to split the collision kernel as 
in [241 Section 2.4]. For the truncated and regularized part Qg (we use the notation 
introduced in step 2 the proof of Proposition 12. ip . [24^ Corollary 2.3] together with some 
basic interpolation yield the following control: 

/ (Q+(S,D)+Q+(D,S))Ddv < Cp 1+2e \\D\\ 2 - 2e 
Jr n 

for some explicit C > and 9 £ (0,1). For the remaining term Q^, we use the same 
control as in [24, Proof of Proposition 2.5] to get 

/ (Q+(S, D) + Q+(D, S)) Ddv<e Ud\\l\ + PIU?J PIU? /a 
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for some e which can be taken as small as wanted by the truncation. Gathering these 
estimates, we get 

Ve>0 f Q + (S,D)Ddv <e\\D\\ 2 T2 + C e 

where C £ depends on weighted L 1 and L 2 norms of S, on L 1 norms on D and on e. Using 
equation (I2,30p with i = 0, Lemma 12.31 to treat the term L(D), and some elementary 
interpolation, we deduce that 

||D||r2 <C||D||rl 

for some constant C > 0, which concludes the proof for k = using the previous step on 
the L 1 moments. 

Step 3. Control of the H k norms. From the previous step and some interpolation, in 
order to conclude it is enough to prove (|2.26p for any k £ N and m = 1. We proceed by 
induction on k. For any i £ ~N N , the equation satisfied by d % D is 

d i Q + (S,D) + d i Q + (D,S) -d l (L(D) S) - L{S) d i D 



~ Z) ( i )d t 'L(S)d t ~ i 'D-p(l-a)d t V-(vD) = 0. 

0<i'<i ^ ' 

We deduce that 

(2.30) cf (&D) 2 (l + \v\)dv< [ (d i Q + (S,D)+d t Q + (D,S))d i Ddv 



d i 'L(D)d i - i 'Sd i Ddv 

&'L(S)dHD&Ddv 



0<i'<i 



0<i'<i 



droping the non-positive term. 

The induction is initialized by Step 2. Let us assume the induction step k > to be 
proved, and let us consider some i £ N w such that |i| = k + 1. Using equation (12.301) 
and |24} Theorem 2.5] to estimate the gain term, we find easily 



\\d l D\\ L 2 < C [\\D\\ L x + \\D\\ H , +(3 - N)/ i 

for some q > 0. Therefore we obtain by interpolation (since (3 — N)/2 < 1 for N > 2), for 
another q' possibly larger: 

\\D\\ H k+i < C (\\D\\ L i + \\D\\ Hk 

This concludes the proof, using interpolation, the induction hypothesis k, and the Step 1 
on the L 1 moments. □ 
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3 The elastic limit a — > 1 

3.1 Dependency of the collision operator according to the inelasticity 

In this subsection we show that the collision operator depends continuously on the inelas- 
ticity coefficient a G [0, 1]. Since it is an unbounded operator, this continuous dependency 
is expressed in the norm of the graph of the operator or in some weaker norm. We start 
showing that this dependency of the collision operator is Lipschitz, and even C 1,v for any 
rj G (0, 1), when allowing a loss (in terms of derivatives and weight) in the norm they are 
expressed. Let define the formal derivative of the collision operator according to a by 

Q' a (9,f) :=V„- ( ( I g('v4a))f('v(a))b\u\ f "" 1 " 1 ^ dadv. 

\JrnJ § n-i V 4 « / 

or by duality 

(Q' a (9,f)^)-= f f f 9*fb\u\ ( lula ~ U ) V^(v' a )dadv*dv. 

Proposition 3.1 Let us fix a smooth exponential weight m = exp(— a \v\ s ), a G (0, +00), 
s G (0,1). Then 

(i) For any k,q G N the exists C G (0, 00) such that for any smooth functions f,g (say 
in S(WL ) ) and any a G [0, 1] there holds 

(3.1) HQa(0»/)llwM (m -i) ^ °k,m ll/ll w^^-l) W^W^m-l) 

(3-2) WJj3> f)\\w*>\ m -i) - Ck > m ll/llw*^ 1 ' 1 ^- 1 ) ^w^im-iy 

(ii) Moreover, for any smooth functions f,g and for any a, a' G [0, 1], there holds 

WQUg, f) - Q%is, /)-(«- oi) Q' a (g, f)\\ w-^^ 

(3-3) <\a- a'\ 2 ||/|| L i +3 ( m -i) ||fflLi +3 ( m -i)- 

(Hi) As a consequence, there holds 

(3-4) \\Q+,(g,f) - Qt(g, f)\\w*(m-i) <C\a-a'\ ll/ll^fc+3,1, „ \\9\\ w 2k+3,i (m -i V 
and for any r] G (1,2), there exists k v G N, q v G N and G (0, 00) such that 

\\Qt{g, f) - Qt>(9, /)-(«- «') Q'*(9, f)h Hm -i) 

(3-5) < C v \a - a'l" ll/H^V 1 ) ||3||^,i (m _ ir 

Proof of Proposition 13.11 First by classical convolution-like estimates (see for instance [28J 
in the elastic case, and |17] in the inelastic case, as well as the proof of Proposition 
below) we easily have (13. ID and (13.11) . 
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Next, in order to prove (|3,3p we proceed by duality. Let us consider ip G 5(1^) and 
define ip := (p(v) q m . We compute 



[Qi(9, f) - Qt>{9, /)-(«- ol) Q' a ( g , /)] Hv) dv 

R^ 

= / H / U(t/«) - ^KO - (a - a') f M^T) 



Hence, if one denotes by ^^.^(a:) := ip(v' a ) (for given fixed values of v,v*,a), we obtain 
(omitting the subscripts for clarity) 

|/| = / \u\bg*f [£(a) — — (a — a') £'(a>)] dvdv*do~ 

< (a — a') 2 / |u| bg* f sup ^"(a)! dv dv* da. 

Jm^xR^xS^- 1 ae(0,l) 

We then easily conclude that (|3.3|) holds using that (v') q (m/) -1 < C (v) q (m)~ l {v*) q (m*) -1 
for some constant C G (0, oo). 

Last, we prove (|3.4p by using the following interpolation on J = Q~^(g, f) — Q^i(g, f) — 
(a-a')Q' a (g,f); 

H J HwtV-i) ^ l^lliy- 2 - 1 ^-!) II J H wj^'^m-i) 

and using (|3.3[) on the first term in the right-hand side, and f|3. 1I3.2|) on the second term 
in the right-hand side. It yields 

\\Qt(9,f) ~ Qi(9,f) ~{a- Oi)Q' a (9j)\\ w y( m -i) < C\a-a'\ ||/|| W ra*+a.i (m -i^ hWw^t 3 ' 1 ^- 1 ) 

and (|3.4p follows by using (|3.2p again. 

Then the proof of (|3.5p is done in the same way using suitable interpolation. □ 

We next state a mere (Holder) continuity dependency on a, which is however stronger 
than Proposition 13.11 in some sense, since it is written in the norm of the graph of the 
operator for one the argument. 

Proposition 3.2 For any a, a' G (0,1], and any g G Lj;(m _1 ), / G W^' 1 ^' 1 ), there 
holds 



(3.6) 



\\Qt(g, f) - Qt>(g, /)ILi( m -i) < e(a - «') ll/ll^'V- 1 ) Il^l^m-i)' 

\\Qa(f,9) ~ <9i(/»5)ILi( m -i) < e(a - II/II w m (to -i) H^ILjCm- 1 )- 

1 

where e(r) = Cr 3+4 / s /or some constant C (depending only on b). 



Proof of Proposition \3.2l For any given v,v* G R , u> = w + 7^ and cr G S 1 ^ 1 
we define x G [0, 7r/2], cos x := | cr • Let us fix (5 G (0,1), R G (l,oo) and let 
us define 9 5 G W 1,0O (-l,l) such that 9 s (s) = 1 on (-1 + 25,1 - 25), 6 s (s) = on 
(-1 + 5,1 - <5) c , < 6 S < 1, |0£(s)| < 3/(5, G R (u) = @(\u\/R) with 6(x) = 1 on [0,1], 
Q(x) = 1 - x for x G [1,2] and 0(a?) = on [2, 00), A(<5) := {cr G S N ~ l ; sin 2 x > <5}, 
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B(S) := {a G S 1 *- 1 ; cos6> G (-1 + 25, 1 - 25) c or sin 2 X < We then split Q + in three 
terms, namely 

Qa Qa Qa' Qa' 

where Qa' r is defined by (j 1 .3[) with b replaced by b r := b6s(a ■ u) @r(u), where Q a ' v is 
defined by (|1.3p with b replaced by b v := b 1a(5) (1 ~~ ®r(u)) and where Qa' a is defined by 
KE3D with b replaced by b a := b (1 - 9 s (a ■ u)) @r(u) + b (1 - 8.r(u)) 1^(<5)- We split the 
proof into three steps. 

Step 1. Treatment of small angles. There exists a constant C G (0, oo) such that for any 
a € (0, 1] and 5 G (0, 1) there holds 

llQi' a (^¥')ll L i( ni - 1 )<C5||V'lliJ(m- 1 )llyll^(^- 1 )- 
Indeed let us consider some £ G and let us proceed by duality. We estimate 

< Il6 a l 



and we conclude using that 1 1 b a 1 1 (L 1 ^- 1 )) < C (5 + max,,^, |i?(<5)|) < C 5. 

Step 2. Treatment of large relative velocities. There exists a constant C = C a>S: t, G (0, oo) 
such that for any a G (0, 1] and 5 G (0, 1) there holds 

C 

(3-7) \\Qi' V &, vOllntm- 1 ) ^ ^277 ll^lliKm- 1 ) llylL^m- 1 ) 

We need the following lemma, which we state below and prove at the end of the subsection. 

Lemma 3.3 For any 5 > and a G (0, 1), there holds 

(3.8) cr€S N ~ 1 , sin 2 x > ^ implies m~ l (v') < m (v ) m~ k (v*), 

with k = (1 - V160) S / 2 . 



In order to prove (j3.Tf) we fix £ G L°° and we argue by duality again. We estimate 
thanks to Lemma I3JJI 



Q^ ,v (tp , (p) £(v) m (v)dv = / \u\ V ip# <p£' (m!) l dvdv*do- 

Jr n xR n xS n ~ 1 

< — I \u\ 2 b v ip*<p£' (m)~ k (m*)~* dvdv* da 
R JlR^xR^xS^- 1 

< -r \\t\\L°° IH-I m l - k {.)\\ 2 L ^ \\ip\\ L i im -x)\\<p\\L{( m -i), 



R 

from which we easily conclude since x 1— > xm l ~ k (x) is uniformly bounded by C a ^ s (1 
k)~ l l s , C aiS G (0,oo). 
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Step 3. The truncated operator. Let us prove that there exists a constant C G (0, oo) 
such that for any 5 G (0, 1), a, a' G (0, 1] and R G (1, oo) there holds 

(R? R \ 
\\9\\v-(m-i)\\f\\wi>Hm-iy 

We closely follow the proof of [23j Proposition 4.3]. We consider some £ G L°°, f,g G 
), we proceed by duality and next conclude thanks to a density argument. We have 



I :-. 



[Qt l ' r (9,f)-Qt; r (g,f)]m- 1 £dv 

\u\Q R (u)bs g* f [((v' a ) m" 1 (v' a ) - £(v' a ,) vrT 1 (v' a ,)] dvdv^da. 



With the notations of Lemma [2. 2 1 we perform the changes of variables v ^ v' a = 4> a (v) and 
v t— ► v' a , = <fi a '(v) (f° r fixed v* and a) with jacobians J a and J a >. Observing that without 
restriction we may assume a < a' and therefore O a = + £l Ua (s) C O q ' = v* + f2 w ,(5) 
since s 1— > w s (0) is an increasing function, we get 



I 



+ 



+ 



g* f (m" 1 )' F(0 J" 1 dv' dv* da 
'xS^- 1 Jo a \O a , 

g* £' (m- 1 )' F^- 1 ) - J' 1 ] dv' dv, da 



'xS N ~ 1 JO 



g* £' (m x )' 



J1 + I2 + /3, 



J 1 di/ dv* da 



with F(w;) := \w — w*| 0_r(w — u*) /(w) £>s(<7 • io — u*). For the first term I\ we use the 
backward change of variables v' 1— ► i> = 0~ x (i/) (for fixed u* and <r) and we get 



I1 



'xS N ~ 1 Jr n 



\u\ Qr(u) f g„ £(v' a ) m (v' a ) bs lo<u-a< v dv* dv da 



with 7] := w a ' o u! a r(5) < C <5 -3 / 2 \a — a'| for some constant C G (0, 00). Since u 1— ► |u | S//2 
is an increasing subadditive function, we also have |t>^,| s < (\v\ 2 + |v*| 2 ) S//2 < \v\ s + |v*| s , 
which implies m(v' a ) < Crn. -1 ?™^ 1 for some constant C G (0, 00) (depending of (). As a 
consequence, we obtain 

III I < CR5~ 3/2 



a - a I ||0||£oo ||£|| L oc ||/|| L i (m -i) ||5llLi( m -i)- 

For the term i2> using the backward change of variable v' 1— > v = <p~, (v') (for some 
fixed v* and a) and using the bounds (|2.4p on J a and | J^ 1 — J"/], we obtain 

|7i| < CR5- 3 \a - a'\ \\b\\ L ~ \\£\\ L ~ ||/|Ui (m -i) Wgh^y 

In order to estimate ^3, we introduce at ■= (1 — t) a + t a 1 and, thanks to (|2.3p - (|2.2p . 

we get 



T\<r C \ 'I 
< — \a — a \ 



Ju N xS N - 1 JO a 



\g,\ \£'\{m- l )'\v' -v\ V w F(^(v>)) 



dv 1 dv* dadt. 
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Using finally the backward change of variable v 1 
(12. 4|) on J at , t G [0, 1], on v* + we get 



v ) and the uniform bound 



f R 2 R\ 

\h\ < C ( + J |a-«'l ll&lliy 1 .- IbllLHm-i) 



Gathering the estimates established in Steps 1, 2 and 3, we deduce the first inequality in 
(|3.6p . The second inequality in (|3.6p is proved in a similar way (using symmetric changes 
of variable, allowed by the truncation). □ 

Proof of Lemma 13.31 We proceed in three steps. 

Step 1. Assume first that (2/\/5) \v*\ < \v\ < (\/5/2) Using the fact that x 1— ► re 5 / 2 
is an increasing and subadditive function, there holds 

\v'\ s < (M 2 + H 2 ) s/2 < (9/4) s / 2 \v*\ s , 

and then by symmetry and because s < 1 

i«T < ^ ( 9 / 4 ) s/2 (M s + \ v *\ s ) < I (M s + \ v *\ s )- 

In that case, (|3.8j) holds with k = 3/4. 

Step 2. We shall first show that for any v,v* G and a € S 1 ^^ 1 , there holds 



(3.9) I 
We recall the formula 



.'1 2 |/|2 



f I J6 j"* Z I 



1 + Q 



sin x v + w * 



. v + v* 1 

w := — 

2 2 



1 — a 1 + a 
u H u °" 



v + v* 1 
2 2 



1 - a 



1 + a 



U <7 



Straightforward computations yield (denoting S" = v + -u*) 
'l + a 2 



k/| 2 <^ + 



2 I 
! ' 4 



2 1 — a 2 2 
kt H u cos ( 



H — (5 •«) H — |5| |u| cosx- 



We deduce the bound from above 



\v'\ 2 < L j- + L j r + \S\ \u\ + \S\ \u\ cosx- 



Then by applying twice Young's inequality 



/|2 < ISfiU 1 — 01 ^ ' 



1 1 — a\ 1 + a . . . 
+ ( ^ + — ^— J + — ^— \S\ \u\ cosx, 



< ^(4 + 



2 / 1 , 1 ~ a \ , 1 |2 



1 1 — a 1 + a 
+ ""- , 4 + — + — 



, 1 + a I C|2 2 



< 



2 2! 



|5| 2 (cos 2 X -l), 



from which we deduce (13.91). 
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Step 3. Assume that sin 2 \ — $ an d that either (2/\/5) \v*\ > \v\ or \v\ > (y/b/2) \v*\. In 
the first case, we have 

\v + v*\ > (1 - (2/V5)) k| + (2/V5) H - M > (1 - (2/VE)) 
which then implies 

+ > (1- (2/^5)) (V5/2) |v| > (1 - (2/\/5)) |u|. 
The same inequalities are proved in a similar way in the second case. We deduce 

\v + v*\ 2 >^(l-(2/V5))(\v\ 2 + \v*\ 2 ). 

We then deduce from flSJlJ) that \v'\ 2 < (1-5/160) (M 2 + H 2 ) and we conclude that (pTHj) 
holds as in Step 1. □ 

3.2 Quantification of the elastic limit a — > 1 

We begin with a simple consequence of Proposition 13. 11 

Corollary 3.4 There exists ko,qo £ N sue/t that for any a% £ (0, oo) % = 1, 2, 3, there 
exists an explicit constant C £ (0, oo) such that for any function g satisfying 



|£tf,a(0)-£>tf,l(0)l <C7(l-a), 
where we recall that Dn, a is defined in $2.18\) . 

Proof of Corollarv \3.4l We write 

DH,a{g) ~ D H>1 (g) = j j j h \ u \ 9a 9* a ~ 9 9*] dv dv* da (=: h) 

+ j j j b\u\gg*[ log g' a + log g' #a - log g - log dv dv* da (=: J 2 ). 

For the first term, thanks to Proposition 13.11 we have 

|/i| < \\QZ(g,g)-Qt(g,g)\\ L i<C(l-a)\\g\\ 2 w s, 1 . 
For the second term, we write 

\h\ = 2\((Q+(g,g)-Q+(g,g))(v) s ,(v)- 8 logg) 

< 2\\Q+(g,g) - Qi(g,g)\\ Ll \\(v)- 8 logg\\ L ~ 

< C(l- a) \\g\\ 2 w ^ (I log y|L«| + |]oga2|+03) 

< C (1 — a) a\ (| log oi| + | log a?\ + Ca 3 ), 
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thanks to Proposition 13.11 and the bounded embedding H k ° n L l qQ C L°° n W x { for ko,qo 
large enough (see Proposition IB. lh . We conclude the proof gathering these two estimates. 
□ 

Let us recall now two famous inequalities, namely the Csiszar-Kullback-Pinsker in- 
equality (see OE2]) and the so-called entropy-entropy production inequalities (the ver- 
sion we present here is established in [31]) that we will use several time in the sequel. 

Theorem 3.5 (i) For a given function g 6 h\, let us denote by M[g] the Maxwellian 
function with the same mass, momentum and temperature as g. For any < g G 
L\(R N ), there holds 

(3-10) \\g-M[g]\\ 2 L1 <2p(g) [ g\n-^-dv. 

ii" M [91 

(ii) For any e > there exists k £ , q £ € N and for any A £ (0, oo) there exists C £ = 
Ce.A £ (0, oo) such that for any g £ H ke n such that 



g{v)>A- x e- A \ v \ , N|**n24<4 

there holds 

(3.11) CM9) 1 -' {j^^Y' <D«M- 



We have then the following estimate on the distance between G a and G\ for any 
self-similar profile G a . 

Proposition 3.6 For any e > there exists C £ (independent of the mass p) such that 
(3.12) Va€[a ,l) sup ||G a - Gi||jr,i < C e p (1 - a)& 

where we recall that G\ is the Maxwellian function defined by U.25\) ~ (1.27 ). 



Proof of Proposition 13.61 On the one hand, for any inelasticity coefficient a € [ao, 1) and 
profile G a , there holds from (|2.19p together with Corollary 13.41 and the uniform estimates 
of Proposition 12.11 

(3.13) D H ,i(G a ) < D H , a (G a ) + p 2 Oil -a)<p 2 0(1 - a). 

On the other hand, introducing the Maxwellian function Mg with the same mass, momen- 
tum and temperature as G a , that is Mq given by (jl.26p with u = and 6 = 8{G a )/p, 
and gathering (|3.13p . (|3.1ip . (I3.10p with the uniform estimates of Proposition 12.11 and 
interpolation inequality, we obtain that for any q, e > there exists Cq )£ such that 

(3.14) Vq G [a , 1) \\G a - M e \\ 2 + e < C q , £ p 2+£ (1 - a). 

i 

Next, from (|2.16j) . we have 

b\ f f G a G a *\u\ 3 dvdv* — p f G a \v\ 2 dv = (1 — a) f f G a G a *\u\ 3 dv dv* 

Jr n Jr n Jr n 2 J^nJ r n 
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and then 

(3.15) <d\\G a - Me\\ L i+C 2 p 2 (l- a), 

where we have used that G a and Mg are bounded thanks to Proposition 12.11 and we have 
defined 

(3.16) f(0) = p I M e \v\ 2 dv-h / M e M e * |«| 3 dv dv*. 
By elementary changes of variables, this formula simplifies into 

¥(0) = kx e - k 2 3/2 

with k\ = p 2 N and, using (|A.3|) . 

k 2 = p 2 h [ Mx, ,x(Mx, Q ,x)*\u\ 3 dvdv* = 2 3 / 2 p 2 bxm 3/2 (M lfijl ). 

JR N xR N 

We next observe that E C°°(0,oo) and is strictly concave. It is also obvious that the 
equation ^(0) = for > has a unique solution which is 6\ defined in ([1.270 . and that 
we have 

*(0) < &(pi) (0 - 0i) = -k x (0 - 0i)/2 

as well as 

(3.17) * (0) =e[kx- k 2 e 1/2 ] = k 2 e [o\ /2 - e 1 ' 2 ]. 

Plugging this expression for into (|3. 15 j) and using the lower bound ([2.170 on the tem- 
perature and the estimate ()3,14p we obtain that for any e > there is C £ E (0, oo) such 
that 

(3.18) Vae(a 0) l) 1/2 - e\ /2 ^ < C £ (1 - a). 

Namely, we have thus proved that the temperature of G a converge (with rate) to the 
expected temperature 0\. In order to come back to the norm of G a — Gx, we first write, 
using Cauchy-Schwarz's inequality, 

||G„-Gi||^ < \\G a -M e \\ Ll _ N + \\Me-Gx\\ L i_ N 

(3.19) < \\G a -M e \\ L i+C N \\M e -Gx\\ L 2, 

and we remark that 

(3.20) \\M g - Gi|| 2 2 < Cp 2 (0 1 / 2 - 9\ /2 \. 

Gathering (pTT9|) with (pT20|) . (pUg]) and (pTEl]) we deduce that for any e > there is 
C £ E (0, oo) such that 

Va€(a ,l) \\G a - GxW 2 ^ < C £ p 2+e (1 - a), 

— N 

and (|3.12|) follows by interpolation again. □ 



36 



4 Uniqueness and continuity of the path of self-similar pro- 
files 



4.1 The proof of uniqueness 

Theorem 4.1 There exists a constructive ct\ G (0, 1) such that the solution G a of U.cl(J\) 
is unique for any a G [a±, 1]. We denote by G a this unique self-similar profile. 

That is an immediate consequence of the following result. 

Proposition 4.2 There is a constructive constant n G (0, 1) such that 

G,HeG a , ae (1-77,1) 

\\G-G 1 \\ Ll < v , \\H-G 1 \\ Ll <r l 



implies G = H. 



Proof of Theorem 14.11 Let us assume that Proposition 14.21 holds. Then Proposition [37 
implies that there is some explicit e G (0, 1) such that for a G (1 — e, 1] one has 

sup \\G a - Gi\\ L i < T) 

where n is defined in the statement of Proposition 14.21 Up to reducing rj, it is always 
possible to take rj < e, and the proof is completed by applying Proposition 14.21 □ 
Proof of Proposition \4.2\ Let us consider any exponential weight function m with 
s G (0,1), a G (0, +oo), or with s = 1 and a G (0, do) small enough. Let us also define 
C = Co,o,o H L 1 (m^ 1 ) the subvector space of of L 1 (m _1 ) of functions with zero energy, 
ip = C {\v\ 2 — N) Mi ; o,i such that £(ip) = 1, and LT the following projection 

n : LV" 1 ) -»■ O, n(g)=g-£(g)ip. 

Finally, let us introduce $ the following non-linear functional operator 

$ : [0, 1) x (Wj'^m -1 ) n C P)0 ) IxO, 

and 

$(V) : {L\{m~ l )r\C pfl ) -^IxO, 

by setting 

$(«,#) = (il + a)D s (g)-2p£{g), n Q a (g,g) -r a div,,(w 5 ) ) . 

It is straightforward that $(a,G a ) = for any a G [ojo, 1] and G Q G Q a , and that the 
equation 

$(1,9) = (0,0) 

has a unique solution, given by g = G\ = M ^ defined in fjl .25|) . (jl.27p . 

The function $ is linear and quadratic in its second argument by inspection, and easy 
computations yield the following formal differential according to the second argument at 
the point (1, Gi): 

(4.1) D 2 <5>{l,G 1 )h = Ah:= ( AD s (G 1: h) - 2p£{h), 2Q 1 (G 1 ,h)\ 
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where Q a is defined in (JT78J) and 

D £ (g,h) :=h I f 

J JR 



g h* \u\ dvdv*. 



Notice that we can remove the projection on the last argument in (14. ip since the elastic 
collision operator always has zero energy. 

Then we have the 



Lemma 4.3 The linear functional 



A : Li{m- r ) 
h 



1x0 

Ah = L> 2 $(l,Gi)/i 



is invertible: it is bijective with A 1 bounded with explicit estimate. 

Proof of Lemma \4.3[ Since the spectrum of the linear operator L\ defined on L 1 (m _1 ) 
(with domain L\{m~ 1 )) includes as a discrete eigenvalue associated with the eigenspace 
Ker£i = SpanjGi, V\ G\, . . . , vn G\, \v \ 2 G\\ by [27J, Theorem 1.3] and since moreover 
O D Ker£i = {0}, we deduce that it is invertible from O n L](m _1 ) onto O. Moreover 
the work [27l Section 4] provides explicit estimates on the norm of its inverse. We deduce 
immediately that C^ 1 maps O onto itself with explicit bound. 
For any h E L 1 (m _1 ), we decompose 



h = hi 4>\ + h , with h\ :- 



€ K, h E O, 



where we recall that <f>i is defined in (jl.33p . Then, using the characterization (jl.30p of G\ 



Ah=\bi I Gi(v)h (v*)\u\ dvdv* 

lR N xR N 



+hi 



b\ / \v\ 2 Gi Gu \u\ 3 dv dv* — 2p / Gi\v\ dv 

R N xR N JR N 



The claimed invertibility follows from the fact that C* = 2N p 2 8 2 ^ 0. Indeed, from 
(fOj) and (TO|> there holds 



C* := bi / \v\ 2 GiG u \u\ A dvdv* - 2p / Gi\v\*dv 

JR n xR n 



n 2 n 2 

p y l 



n 2 ft 2 
P ^1 



jl/2 



biB^ / Mi ) o,i(Mi )0 ,i)*|ur|u| d dudu 1| ,-2 / M i,o,i M dv 

R N xR N 



bi 0\ 12 v / 2(2iY + 3)m 3/2 (M li0 ,i) - 2N{N + 2 



and we conclude thanks to formula (11.271) . 



□ 
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Let us come back to the proof of Proposition 14.21 We write 

G a -H a = A' 1 \a G a - <3?(a, G a ) + 3>(a, H a ) - AH D 

(4.2) = A- l (hJ 2 ) 

with (recall that the bilinear operators and Q a are symmetric) 

f /i :=4D £ (G 1 ,G a -H a )-(l + a)D(G a ) + (l + a)D(H a ) 

\ i 2 :=n/ 2il + n/ 2i2 

and 

'■= 2Qi{Gi,G a — H a ) — Q a (G a , G a ) + Q a (H a , H a ) 
h,2 := P {l-a)V v -{v (H a -G a )). 

On the one hand, 

h = 2 D(2G 1 - (G a + H a ),G a - H a ) + (1 - a) D(G a + H a , G a - H a ) 

so that 

\h\ < C 3 (||Gi-G a || L i + ||Gi-fla|| L i 

+(1 - a) \\G a \\ L i + (1 - a) \\H a \\ L i) \\G a - H a \\ Ll 

(4.3) < rjti^WGa-H^^ 

with 771(a) — ► when a — ► 1 (with explicit rate, for instance 771(a) = C\ (1 — a) 1//3 ) because 
of Propositions 12.11 and 13.61 
On the other hand, 

^2,1 = Qi{G\,G a — H a ) — Q a (Gi,G a — H a ) + Qi(G a — H a ,Gi) — Q a (G a — H a ,Gi) 

+Q a {G\ — G a , G a — H a ) + Q a (G a — H a , G± — H a ). 

From Proposition 13.21 there holds 

\\Qi(Gi,G a - H a ) - Q a (G\,G a - H a )\\ L i( m -i) < e(a) \\G a - H a \\ L i^ m -i^ 

\\Qi\G a — H a , Gi) — Q a {G a — H a , Gi)||ii( m -i) < e(a) \\G a — -#a||£i( m -i) 

with e(a) — > as a — > 1 (with again explicit rate, for instance e(a) = C[ (1 — a) 1 / 12 if 
s = 1/2 in the formula of m). From elementary estimates in L 1 (m _1 ) we have 

\\Qa{Gi — G a , G a — H a ) + Q a (G a — H a , G\ — H a )\\ L i^ m -i^ 

< C 4 (\\G a - Gil^i^-i) + \\H a - Gi|| L i( m -i)) \\G a - Ha|Li( m -i)- 

Together with Propositions 13.61 we thus obtain 

(4-4) \\h,l ||ii( m -i) < 772(a) \\G a - i?a|| L i (m -i) 
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for some 772(a) — ► as a — > 1. Here we can take for instance (when s = 1/2 in the 
formula of m) 7? 2 (a) = C 2 (1 - a) 1/12 for some C2 6 (0, 00) by picking a suitable e and 
interpolating. 

Finally from Proposition 12.71 there holds 

(4.5) Il^2,2||z,i( m -i) < C$ (1 - a) ||Gq, - F Q || L i (m -i). 
Gathering (|4.3p . (|4T4]| and (|43|) we obtain from (|4T2j) and Lemma liTBI 

||Gq, - H a ||ii( m -i) < Vi a ) \\A X \\ \\G a - #a||zi( m -l) 

for some function r/ such that 77(a) — > as a — > 1 (with explicit rate). Hence choosing a\ 
close enough to 1 we have 77(a) < 1/2 for any a G [ai, 1). This implies G a = H a 

and concludes the proof. □ 

4.2 Differentiability of the map a 1— > G a at a = 1 

Lemma 4.4 T/ie map [ai, 1] — >• L 1 (m _1 ) ; a 1— > G a is continuous on [a\, 1] and differen- 
tiable at a = 1. More precisely, there exists G' x E L 1 (m _1 ) and for any rj £ (1,2) there 
exists a constructive E (0, 00) such that 

(4.6) ||G a -Gi-(l-a)Gi|| £ i (m -i ) <C r „(l-a)'' Vae(a ,l). 



Proof of Lemma 14.41 We split the proof into four steps. 

5*tep 1. For the continuity we use a classical stability argument. Let us consider a sequence 
(a„)„>o such that a n E [ai, 1] and a n — » a. From the uniform bound (12. ip . we may extract 
a subsequence (G Qii ,) which strongly converges in L 1 (m~ 1 ) to a function G a . Passing to 
the limit in the equations (11.301) associated to the normal restitution coefficient ot n and 
written for G« n ,, we deduce that G a satisfies f j 1 . 30 j) associated to the normal restitution 
coefficient a. From the uniqueness of the solution proved in Theorem 14.1^ there holds 
G a = G a and thus the whole sequence G an converges to G a . 

Step 2. We next prove that there exists an explicit constant C such that 
Va G [a<i,l] \\G a - Gi\\ L i( m -i) < C(l - a). 

We write 

G a — G\ = A' 1 [AG a -^{a,G a ) + ^{l,G 1 )- AG X ] 
(4.7) = A^iJ^h) 

with 

Ji := 4De{Gx, G a - G x ) + 2D e (G 1 , G 1 ) - (1 + a) D £ (G a , G a ) 

J2 :=nj2,i+nj2,2 
and 

^2,1 := Qi{Gi,G a ) + Qi(G a , G\) — Q a (G a , G a ) 
J2,2 :=p(l - a) V„ ■ (v(G a )). 
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On the one hand, 

J x = -2D e (Gi - G a , Gi - G a ) + (1 - a) D(G a , G a ) 

so that 

\Ji\ < C\\Gi - G a \\ 2 Tl +(7(1 -a). 

On the other hand, 

J2,i = — Qi(G\ — G a , G\ — G a ) + Qi(G a , G a ) — Q a {G a , G a ). 
Hence using Propositions 12.71 E5H and the bound (|2.ip . we deduce 

\J2,i\ < ||G a -Gi||^ + C(l-a) 

and we also have straightforwardly J2,2 = 0(1 — a). Gathering all these estimates, we 
thus obtain from (|4.7p 



\G a - Gi[| L i( m -i) < P *| 



|G Q - Ci lli^m- 1 ) + C(l - a) 



Using then the explicit result of quantification of the elastic limit in Proposition 13.61 we 
have that for some a 2 G 1) close enough to 1: 

Vae [a 2 ,l] WA^WWGa-GiW^^ < - 

and thus we get 

VaG [a 2 ,l], \\G ol -G x \\ L i {m - i) < 2(7 (1 - a) 

which implies the claimed estimate. 

Step 3. In order to prove the differentiability we must slightly improve the estimate estab- 
lished in the preceding step. On the one hand we exhibit what should be the derivative 
of G a at a = 1, and denote it by R. Formally differentiating equation (|1.30p at a = 1 we 
have 

Q' 1 (G 1 ,G 1 )+2Qi(R, Gi) + pV v ■ (vG x ) = 0. 
On the other hand, we may compute 

(<X(Gi,Gi),|.| 2 ) = - III b\u\G 1 G u (\u\a-u)- (\u\a)dvdv*da 

4 J J J R N xR N xS N-l V / 

(4.8) = 2D £ {G 1 ). 

Next, diving the equation (|1.23p on the energy of G a by (1 — a) and formally differentiating 
the resulting expression we get 

2p£(R)-De(G 1 ,G 1 )-4D e (R,G 1 ) = 0. 

We now rigorously define R in the following way 

G' x = R := A _1 ( - De(Gi,Gi),—F), F := Q' a (Gi, G x ) + p V v ■ (v G x ). 
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Note that R is well-defined since £{F) = because of f|4.8[) and the definition of G\. 

Step 4. We finally come back to the Step 2 and we shall construct a Taylor expansion of 
order 1. We want to estimate 

G a - Gx + (a - 1) G[ = A' 1 (j x - (a - 1) D £ (G U Gi), J 2 - (1 - a) f) . 
On the one hand 

Ji - (a - 1) Dg(Gi,G{) = -2 De(G\ - G a , Gi - G Q ) + (1 - a) (D(G a , G«) - D £ (Gi, Gi)) 
so that we obtain straightforwardly 

|Ji-(a-l)A:(Gi,Gi)| <G(l-a) 2 . 

On the other hand, 

J 2 - (1 - a) F := n J 2j i + n J 2 , 2 

with 

^2,1 = — Qi(Gi — Gq,, Gi — Gq,) + Qi(G a — Gi, Gq,) — Q a [G a — G%, G a ) 

+Qi(G l ,G a - Gi) - Q a (Gx,G a - Gi) + (1 - a) V, • (« (G Q - Gi)) 

and 

J 2 ,2 = Qi(G 1 ,G 1 ) - Q a (Gi, GO - (1 - a) K 

It is clear from Propositions 13.14 the bound of Step 2, and some interpolation with the 
uniform bounds (|2.ip . that 

ll^2,l||Li(m-i), 11^2,2 lU^m- 1 ) < Cfe (1 ~ a) fe 
for any k G (1, 2). □ 

5 Study of the spectrum and semigroup of the linearized 
problem 

In this section we shall obtain the geometry of the spectrum of the linearized rescaled 
inelastic collision operator for a small inelasticity, as well as estimates on its resolvent 
and on the associated linear semigroup. This is based on the properties of the elastic 
linearized operator and some perturbation arguments again. In order to do so, one needs 
some common functional "ground" for the the linearized operators in the limit of vanishing 
inelasticity. This common functional setting is given by the study [27| in which the spectral 
study of the elastic linearized operator is made in L 1 spaces with exponential weights e a , 
a £ (0,+cxd), s £ (0,1). 

We thus consider the operator 

9 Q a (g,g) -r a V„ ■ (vg) 

and some fluctuations h around the self-similar profile G a : that means g = G a + h with 
h E L 1 (m _1 ) where m is a fixed smooth exponential weight function, as defined in (jl.28p . 
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The corresponding linearized unbounded operator £ a acting on L l (m~ l ) with domain 
dom(£ a ) = W t ' (m~ x ) if a ^ 1 and dom(£i) = L\{m~ l ), is defined in (|1.31j) (it is 
straightforward to check that it is closed in this space). Since the equation in self-similar 
variables preserves mass and the zero momentum, the correct spectral study of £ a requires 
to restrict this operator to zero mean and centered distributions (which are preserved as 
well), that means to work in L 1 (m _1 ). When restricted to this space, the operator £ a is 
denoted by £ a . We denote by R(£ a ) the resolvent set of £ a , and by TZ a (£.) = (A* — 
its resolvent operator for any £ G R(£ a ). 

Let us recall that for the linearized elastic hard spheres Boltzmann equation the spec- 
trum and the asymptotic stability have been studied by many authors since the pioneering 
works by Hilbert [20], Carleman [12] and Grad [H], and we refer for instance to |27j for 
more references. The result established for C\ (and translated straightforwardly to C\) 
in [27] is the following: 

Theorem 5.1 (i) There exists a decreasing sequence of real discrete eigenvalues (fJ, n )n>i 
(that is: eigenvalues isolated and with finite multiplicity) of L\, with "energy" eigen- 
value fii = of multiplicity 1 and "energy" eigenvector (pi (defined in il.33]) ), ^2 < 
and lim/^ n = fi^ G (— 00, 0) such that the spectrum £(£1) of £1 in L 1 (m _1 ) writes 

£(£].) = (—OO, Moo] U {/U n } neN . 

In particular, £1 is onto from O n L}(m~ 1 ) onto O. 

(ii) The resolvent 7£i(£) has a sectorial property for the spectrum substracted from the 
"energy" eigenvalue, namely there is a constructive fj,2 < A < such that 



A=<£€C, arg(£ + A) G 



Willi _ 

— — ,— — and 3f?e£ < — 
4 ' 4 J s ~ 2 

(Hi) The linear semigroup Si(t) associated to £1 in L 1 (m^ 1 ) writes 

Vt>0 Si(t) = Ui + Ri(t), 

where IT is the projection on the eigenspace associated to /ii and Ri(t) is a semigroup 
which satisfies 

Vt>0 ||iM*)lki(m-i)<Ce^ 
with explicit constant C . 



The main result proved in this section is a perturbation result which extends Theo- 
rem 15.11 in the following way. Let us define for any x G R the half-plane A x by 

A, = {£GC, 3fie£>x}. 

Theorem 5.2 Let us fix p, G (/j,2, 0), k, q G N and m a smooth weight exponential function 
with s G (0, 1). Then there exists Q2 G (a%, 1) such that for any a G [012, 1] the following 
holds: 
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(i) The spectrum H(C a ) of C a mWg' 1 (m l ) writes 

E(A») = E a U £ a C A£, 

where [i a is a 1-dimensional real eigenvalue which does not depend on the choice of 
the space Wj'^m -1 ) and satisfies ( TT7 



(mJ T/ie resolvent H a (£) in wj'^m x ) is holomorphic on a neighborhood of A^\{/x a } 
and i/iere are explicit constants Ci,C 2 such that 

and 

C 2 

\\K a (fj, + is)\\ w k+iA, m _ lww fc,i/ m _ii < 



, + 1 (m-i)-»Vf*' A (m-i) - 1+| S |' 

fm) T/ie linear semigroup S a (t) associated to C a in Wg ,1 (m~ 1 ) writes 

S a (t) = e^ t U a + R a (t), 

where U a is the projection on the (1-dimensional) eigenspace associated to \x a and 
where R a (t) is a semigroup which satisfies 

(5-1) II^WHw^^m-iHW^^m-i) ^ C k e ^ 

with explicit bounds. 



Remark 5.3 Note that we do not claim that the resolvent 1Z a is sectorial for a < 1 and 
it is likely that indeed it is not (because of the contribution of the drift term). Moreover, it 
is not clear how to make the spectral study in the Hilbert setting L 2 (m~ 1 ) with convenient 
weight function m. In particular, we are not able to prove Proposition \3.2\ in an I? 
framework. In such a situation the spectral study and the obtaining of constructive rate of 
decay on the semigroup become tricky. Let us emphasize also that ( as most of the results 
established in this paper) this result is not an easy consequence of perturbation theory of 
unbounded operator since the elastic limit a —* 1 is strongly bad-behaved (for instance 
neither the relative bound nor the operator gap of ]21^ go to 0) because of the anti-drift 
term. 



5.1 Recalls and improvments of technical tools from [27j 

Proposition 5.4 In the statement of Theorem I5.il one can replace everywhere L l {m~ l ) 
by Wj'^m -1 ), k,q G N. 

Let us first recall the key decomposition of C\ in [27} Section 2] (re-written within the 
notation of this paper): 

Let 1 e denote the usual indicator function of the set E, let : R — > be an even 
C°° function with mass 1 and support included in [—1, 1] and : R — > IR + a radial C°° 
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function with mass 1 and support included in B(0, 1). We define the following mollification 
functions (e > 0): 

( e e (x) = e - 1 e( e - 1 x), (x g m) 

\ & e (x) = e~ N Q(e~ 1 x), (x£R N ). 
Then we consider the decompositions 

C l {g)=C\(g)-C»(g) with C{g) := ug 

where C\ splits between a "gain" part Cf (denoted so because it corresponds to the 
linearization of Q + ) and a convolution part C* (not depending on a) as 

Ci(9) = £f(g) - C{g) with £*(g) := M [g * $], 

(we do not write the subscript 1 when there is no dependency on a). Then for any 5 G (0, 1) 
we set 



£+ s (g)=l s (v) / <Z>(\v-v*\)b 5 (co S 6)lg'Mi+M'gi}dv*da, 

where 

Is = @s* 1 {|-|<<5- 1 }> 

and 

b 5 {z) = (9,52 * 1{_i+252< 2 <1-2«52}) K z )- 

This approximation induces C\ : $ = C^ s — C* — C u . Then the key result is that this 
approximation converges (in the norm of the graph) to the original linearized operator 
C\ as 5 — > 0, first in the small classical linearization space L 2 (G'^ 1 ) (this technical result 
was in fact mostly already included in Grad's results [18]). and second most importantly 
in the larger space L 1 (m~ 1 ). On the basis of this approximation result the spectrum is 
then proved to be the same in both functional spaces, and then the norm of the resolvents 
within these two functional spaces are related by an explicit control. 

Hence the keys elements of the proof which are to be extended are, on the one hand, 
the approximation argument (which has to be extended from an L 1 (m _1 ) setting to an 
W q ,l {rn~ l ) setting), and, on the other hand the explicit control on the resolvent in the 
space L 2 (M _1 ) provided by the self-adjointness structure of the collision operator in this 
space and the explicit estimates on the spectral gap (see [6]), which has to be extended to 
an H k (M~ l ) setting. Then the rest of the proof of [27] would extend as well (up to minor 
technical modifications) to W k,p {m~~ l ). 

Therefore for the first point let us prove the 
Proposition 5.5 For any fc,g£N and g G W ^{wT 1 ), we have 

where e(5) > is an explicit constant going to as 5 goes to 0. 

Proof of Proposition \5.5\ The case k = q = is provided by Proposition 13.21 Then 
higher-order derivatives follows by differentiation, and the incoporation of a polynomial 
weight is trivial. □ 
Concerning the second point let us prove the 
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Proposition 5.6 The spectrum of Ci in L 2 (M 1 ) is the same in any 

k G N. Moreover the control on the resolvent, which was (self-adjoint operator) 

II^i(£)IIl2 ( m-i) < ' 



in the space L 2 (M ), extends into 
with 

' £ G C, ara(£ + A) G 



dist(£,E(£i)) 
Cfe 



dist(£,£(£i))' 



A 



3ir 3tt 

T'T 

/or any A; G N and some explicit constant Ck > 0, 



and 5fte£ < — 



Proof of Proposition 15.61 A quick way to prove the result for instance is the following. 
It is easy to prove by induction on k G N the following estimate on the Dirichlet form: 



|s|<fc \|s|<fc 



for some explicit > and a s > 0, |s| < k, and where II denotes the orthogonal projection 
in L 2 (M~ X ) onto the functions with zero mass, momentum and energy. Therefore we 
deduce on L\ that its semigroup satisfies 

VfcGN, \\e t&1 \\ H H { M-i)<Ck 

and that obviously the same is true on the stable subspace of functions with zero energy. 
Then by interpolation with the rate of decay of the semigroup for functions with zero 
energy in L 2 (M _1 ), we deduce that 

Ve > 0, k G N, lle^nlltffc^-!) < C e , fc e-<'"- e > t 

for some explict C e fc > 0, and where II is the orthogonal projection in L 2 {M~ l ) onto 
functions with zero energy. This implies on the resolvent that for any k G N, 

V£GA ||fti(OII/r* ( M-i) < Ci 

with 



A=UeC, arg(£ + A) G 



37T 37T 
T'T 



and 3fte£ < — 



for some explicit CL > 0. Then the result follows by straightforward interpolation with 
the estimates on the resolvent in L 2 (M _1 ). □ 

Then we can conclude to the following extension of point (ii) of Theorem 15.11 

Proposition 5.7 We have 

h,q 



V£ G A, ll^i(0ll w M (m -i) < a k , g + + A , 



with 



A = £ € C, arg(£ + A) G 



3-7T 3"7T 

T'T 

for any k, q G N and some explicit constant ak,q, bk, q > 0. 



and 5fte£ < — 
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5.2 Decomposition of C a and technical estimates 

We fix once for all some fx G (fj,2, 0) and we split the proof of Theorem 15.21 into four steps, 
detailed in the following four subsections. 
Let us introduce the operator 

P a = d - C a = £+-£+ + r a V v ■ (v •). 

Our first step in this subsection is to estimate the convergence to of the first part of 
this operator in suitable norm. Namely we prove 

Lemma 5.8 (i) For any k,q 6N, there exists C = Ck,q, m such that 

(ii) For any k,q G N, there is a constructive function e : (0, oo) —* (0, oo) satisfying 
e{a) — > as a goes to 1 and such that for any g G W^^m^ 1 ) 

|| (£+-£+) (5)||^,i (m _ 1} < e{a) NI^m (m - 1} - 

(Hi) There exists C G (0, oo) such that for any g G W^' l {m~ l ), we have 
||(£i - C a ) (5)IIli (to -i) < C(l - a) \\9\\ w 3,x {m -iy 

Proof of Lemma \5. 81 The case k = q = is proved in Proposition 13.21 Then higher-order 
derivatives are obtained from the L 1 (m _1 ) estimates by straightforward differentiation, 
and the incoporation of polynomial weights is trivial. □ 
Now let us consider some £ G C and let us define 

As = C+ s - C* 

and 

BaAO = " + t + (£t,6 ~ C t) + Pa 

(let us recall that the approximation Ct s was defined in the beginning of Subsection 15.11 
It yields the decomposition 

C a -i = As-B QtS {0- 

Then we have the 

Lemma 5.9 Let us consider any k,q G N and £ such that 5fte£ > — mini/. Then 

(i) For any 5 > 0, the operator As : L 1 — > W^^^m^ 1 ) is a bounded linear operator 
(more precisely it maps functions of L 1 into C°° functions with compact support). 

(ii) For 5 G [0,5*] and a G [ct2, 1] for some constructive 5* > and 012 G (ai,l) 
(depending on a lower bound on dist(£, v(WL N ))), the operator 

B a ,s : W k q +l>\m- 1 ) - W^(m- 1 ) 

is invertible 
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(Hi) The inverse operator B aj s{C) 1 satisfies for 5 £ [0, 5*] and a S [0J3, 1]; 



l^(m->h^( m -') - dist(SRee,^(R JV )) 

and 

||^(e)" 



v — 1 1 1 ^ C2 



for some explicit constants C\,C2 > depending on k,q,5*,a2 and a lower bound 
on dist(5ReC,z/(M JV )). 

Proof of Lemma E3 For £ € u(R N ) c , it was proved in |27l Proposition 4.1, Theorem 4.2] 
the convergence to of — C±) as 5 — > (which was done in L}(m _1 ) — > L 1 (m _1 ) 

in [27] and is extended in any W^^m -1 ) — > Wg'^m -1 ) by Proposition I5.5p . we deduce 
as in [27] that for 5 small enough (depending on a lower bound on the coercivity norm of 
v + £, that is on a lower bound on dist(£, z/(R ))), we have 

||(£+ 5 - /^ffll^i^) < \ + g\\ w ^ ■ 

It was also proved that As maps functions of L 1 into C°° functions with compact 
support (with explicit estimates). 

Let now consider B a s(£) only in the case k = q = (estimates for higher-order deriva- 
tives and weights are obtained by straightforward differentiation and computations) . From 
Lemma [5.8l we have for a close enough to 1 (depending on a lower bound on dist(£, u(M. ))), 

||(£^ -£+) 5 || L1(m _ 1) < - + 

By considering the semigroup on L 1 (m _1 ) of B a $(£) and computing the evolution of the 
norm in symmetric form using the formula for the differentiation of the complex modulus 
of a function 

Vhh + hVh 
V|/l| = ^H ' 

it is easily seen that 

||e B «- i(5) *5|| L i (m -i) > + &e0 gW^^-i) - -\\(v + MeO g\\ L Hm-i) 

and therefore for a close enough to 1 (depending on a lower bound on dist(£, z/(R JV ))), we 
deduce that 

\W BaA0t 9\\Li(m-i) >\\\(v + Me09\\LHm-i) 
and thus that the operator is invertible with its inverse bounded by 

2 

\\B«A0 llii(m-i) < dist(3?ee,KM JV ))' 

Moreover by computing separately the evolution of the L l (m~ l ) norm in non-symmetric 
form (thus keeping v + £ but creating a term of the form 0(1 — a) times a W^' 1 (m -1 ) 
norm) and the evolution of the W i ' 1 (mT 1 ) norm in symmetric form: it yields easily 

\\e Ba ' Si0t 9\\w^(m-l) > 2 W( V + 0^IUi(m-i) + 2 11^ + V»,0|Ui( ro -i) 

which implies the result, by droping the second term. □ 
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5.3 Geometry of the essential spectrum and estimates on the eigenvalues 

First concerning the geometry of the spectrum, following the same strategy as in [27\ 
Subsection 3.2] we can prove the 

Proposition 5.10 Let us pick any k, q G N and m a smooth exponential weight function 
(as defined in hi. 28]) ). Then for any a G [«2,1], the spectrum of C a in wj'^m -1 ) is 
composed of a part included in Al„ (containing all possible essential spectrum), and a 
remaining part included in A_„ exclusively composed of discrete eigenvalues. 

Proof of Proposition \5.10\ We follow the same method as in the proof of [27, Proposi- 
tion 3.4]. One uses the decomposition 

C-a = A$ — B at g(0), 

the compactness of the first part ^4,5 and the coercivity 

||#a,5(0)|Ui( m -i) > ||f 3|Ui( w -i) - e(6) \\v £f||zi( TO -i) 

of the second part (where e(<5) — > as 5 — > 0). Then one applies Weyl's theorem and show 
that (for any 5 > 0) A_ UQ+e ^ has to be a Fredhom set with indices (0, 0) (except possibly 
for a countable family of points) since [a, +oo) is included in the resolvent set for a big 
enoug. □ 

Second concerning the discrete part of the spectrum, that is the isolated eigenvalues 
with finite multiplicity, following the same strategy as in [271 Proof of Proposition 3.5] we 
can prove the 

Proposition 5.11 Let us fix ft £ (//2>0). Then for any a £ [ct2, lj (where «2 is obtained 
from Lemma \5. 9\ for this choice of fi), for any /x G and (p G W 1 ' satisfying 

C a (4>) = ii<f> 

in L l , we have 



W/M( m -i) < Cfc,; 



J 2 

for any k G N and m = exp(— a \v\ s ), a > 0, s G (0, 1), where the constant Ck, m depends 
on k, m and a lower bound on ft — fi. 

Proof of Proposition \5.11\ Let us sketch the idea of the proof. We use the decomposition 

o = c a 4> - ft 4> = A & 4> - B a j(u)(j) 

and the fact that for the choices made for u and a in the assumptions we have (adjusting 
5 as in Lemma l5.9f) B a ^{n) is invertible in any W k,l (m~ l ) with explicit bound, and As 
maps L 1 into C°° functions with compact support. □ 



Remark 5.12 An alternative proof could be to adapt the proof of Proposition 2.7 
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5.4 Estimate on the resolvent and global stability of the spectrum 

Lemma 5.13 Let us pick k, q E N and m a smooth exponential weight function (as defined 
in ll.28\) ) and consider the operator C a in Wj'^m -1 ). Then 

(i) For any £ E R(C\), there is E [02, 1) such that £ E R(C a ) for any a E [ag, 1]. 

(ii) More precisely, the resolvent TZ a {i) satisfies the two following estimates for a E 
[a 2 , 1); 

C\ + C 2 \\T^i(0\\w k +^ 1 {m- 1 ) 



v - 1 - c 3 (1 - a) \\nx(0\\ w ^ {m ^ 



l^-o(C)llw fe + 1 ' 1 fm-l^W*' 1 l'm-M - 



^(m-WV) " 5(0 1 - C 3 (1 - «) ll^lCOIIwJ+M^-i) 

mi/t 5(£) := dist(£, ^(R-^)) and where the constants C^C^, i = 1,2,3 depend on a 
positive lower bound on dist(5fte£, u(M. N )). 

(Hi) Finally, for any compact set K C p{£-i) there exists ax E [a%, 1), Ck E (0, 00) such 
that 

V£ E K, a E (ok, 1] ll^a(Oll w fc.i (m -i) < c k, 

V£ €. K, a, a' E (ajc, 1] ll^aCC) ^ - ^a'(C) ^Ik^m- 1 ) < Ctf (1 _ a ) INIwf' 1 - 
Proof of Lemma \5.13l We split the proof into three steps. 

Step 1. Let us consider the following operator defined from Wg' 1 (m _1 ) to W^l' 1 (m~ l ) 
(which is seen to be well-defined at a glance) 

Some straightforward computations show that 

(£* - l a ,s(0 = -As B^siO' 1 + W + [id - P a Ki(0] A 5 B^siO' 1 
which simplifies into 

First using that 

the control of 7li(£,) in W^j^'^TO^ 1 ) and the regularization property of As we deduce that 

K a ,s(0 = p a fti(0 As BsAO' 1 = 0(1 - a) 

in the norm of bounded operators on W^'^m -1 ), and therefore for (1 — a) small enough 
(with explicit bound) we get that 

WKaAOWw^^-i) < Cs (1 - a) llftiCOIIwJ+J-V- 1 ) K 1 
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and Id — K a ^(£) is invertible in Wq' l {m l ). As a consequence 

(C a - l a ,s(0 (Id - K^Or 1 = Id w *,i (m -i } 

and we have proved that C a — £ admits a right-inverse, namely so that (Id — 

K a> s(£)) . This proves that the operator C a — £ is onto. 

Step 2. In order to show that C a — £ is invertible and that we have identified the 
resolvent it remains to prove that it is one-to-one. Let us consider the eigenvalue equation 
{C a — £) h = which writes 

(A - Oh = Pah 

from which we deduce (using Proposition 15.111 to get regularity bounds on h) 

< C (1 - a) \\K X (0 Hw*'V-i) IHIw^v- 1 ) 

< C" (1 - a) H^lCOIIw^Cm-i) H^llwJ'^m-i)- 

Therefore for (1 — a) small enough (depending on the norm of 7£i(£)) we have that nec- 
essarily h = 0, and thus the operator (C a — £) is one-to-one. 

For a satisfying all the previous conditions, the operator {C a — £) is bijective from 
^fqXx (. m ~ 1 ) to W (? ' 1 (m~ 1 ) and its inverse is given by 

from which we get the desired bound on the resolvent thanks to the study of -B a ,5(C) _1 m 
Lemma 15.91 At this point we have proved points (i), (ii) and the first estimate in (iii). 
Step 3. The second estimate in point (iii) is obtained from the resolvent identity 

Ka{0 ~ ftl(0 = K a (0 [Cx - C a ] Kx(0, 

together with the previous estimates on the resolvent and point (iii) in Lemma 15.81 □ 

Remark that this lemma proves the point (ii) in Theorem 15.21 Moreover, as a con- 
sequence of this estimate on the resolvent H a (£), we may go one step further in the 
localization of the spectrum of C a around 0. 

Corollary 5.14 Let us fix Jx 6 (/X2,0). In any Wj'^m -1 ) there is some constant C G 
(0, oo) such that 

VaG[a 2) l], S(£ Q )n A A c 5(0, £7(1 -a)). 

Proof of Corollary \5.14l The proof follows from the estimates in point (ii) of Lemma [5.131 
together with the fact that (Proposition 4.1 of [27] in L 1 (m~ 1 ) extended to W^'^m -1 ) by 
the previous discussion): 

VeeA^, \\K 1 (0\\ w ^ {m -r ) <a + ^ 

for some explicit constants a, b > 0. We get thus that [|^a(0llw fc,1 (m- 1 ) < 00 if £ ^ ^-p, 
and |£| > C (1 — a), which concludes the proof. □ 
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5.5 Fine study of spectrum close to 

Let us fix r E (0, \p,\] and let us choose any a r E [ot2, 1) such that C (1 — a r ) < r (with the 
notations of Corollary 15. 14p in such a way that n C B(0,r) for any a E [a r , 1]. 

We may then define the spectral projection operator (see [ST]) 

(5.2) n a: =_J_/' ^ a (C)dC 

27 ™ 7s(0,r) 

in any W^'^m -1 ), with 5(0, r) := {£ E C, |£| = r}. The operator H a is the projection 
operator on the sum of eigenspaces associated to eigenvalues lying in the half plane {£ E 
C, 3fte£ > — r}, see 21] • I n particular the operator IIi is the projection on the energy 
eigenline R , where we recall that <j)\ is the energy eigenfunction defined by (jl.33p . 

Lemma 5.15 The operator H a satisfies 

(i) For any k E N and any exponential weight function m (as defined in tl.28\) ). it is 
well-defined and bounded in Wj(m _1 ). 

(ii) Moreover there is a constant C > (depending on m) such that 

(5.3) Va, a' E [a r , 1] \\U a - '^■a>\\ w ^ 1 ( m -i)^ L i( m -i) < c W ~ a \- 

Proof of Lemma \5.15\ It is a straightforward consequence of (15. 2h and Lemma 15.131 □ 

Corollary 5.16 There exists 03 E [02, 1) smc/i i/ia£ for any a E [03, 1) there holds 

Y>(C a ) n Afi = {/x a } and i/ie eigenspace associated to \i a E R is 1- dimensional. 

This eigenvalue is called the energy eigenvalue. We may furthermore remark that Corol- 
lary \5.14\ implies 

(5.4) Vae[a 3 ,l) \/J, a \ < C (1 - a). 

Proof of Corollary \5.16l We already know that £(£<*) n A^ is entirely composed of 
discrete spectrum. Therefore we have to prove that it is of dimension 1. Indeed once this 
is proved, the fact that \x a E R is trivial since the operator is real, and the control (15. 4p is 
trivial from Corollary 15.141 

Let us define the space X a := II a (L 1 (m _1 )) + IIi(L 1 (m _1 )) endowed with the norm 
II ' llL 1 (m- 1 )- From Proposition 15.11] there exists a constant C\ > such that 

VV> E X a , IMIvtf^m-i) ^ ^1 ll^llLi(m-i)- 

Thanks to the definition of Hq, and IT and to Lemma 15.151 we then get 

||TT „ „ ^ n IK^aW-^lW)V'IUi(m-i) 

II Q - IIi hc Q — x a < Q2 sup sup 



< C 2 (l-a) sup — —r 

< C 2 '(l-a)<l, 
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for (1 — a) small enough. By classical operator theory (see for instance the arguments pre- 
sented in [2U Chap 1, paragraph 4.6] in order to prove [211 Lemma 4.10]) one deduces that 
dimension(II Q ) = dimension(IIi). Since dimension(IIi) = 1 (as recalled in Theorem 15. ip . 
this concludes the proof. □ 

Let us introduce for any ip £ L 1 the decomposition 

where niip £ R is the coordinate of Hitp on (defined thanks to the projection IIi). 
For any a £ [«3, 1) we denote by <f> a the unique eigenfunction associated to /J, a such that 
\\4> a \\ L i = 1 and > 0. 

We can now establish a first order approximation of the eigenfunction (f> a . 

Lemma 5.17 For any k, q G N and any exponential weight function m (as defined in \1.28\) ), 

there exists C such that 

(5-5) Vae[a 3 ,l] \\<j> a - ^lll^M^-i) < C (1 - a). 



Remark 5.18 We immediately deduce from Lemma \5.17\ that a (O) < for a close 



enough to 1, and therefore, we get that this definition of <p a coincides with the definition 
in TheoremM. 1[ 



Proof of Lemma \5.17\ On the one hand, from the normalization conditions, we have 



Mi* ~ H n i0a|lLi 



LWaWLl = 1,111 VattLl- 



< 1 1 (pa ~ III <t>n 1 1 r l = I Hi 

We then deduce 

( 5 - 6 ) ||01 - 4>o\\lI < ||ni0 Q - 0a 1 1 2,1 + ||n^ Q || L i < 2 1 1 II] 1 - a ||2,i- 

On the other hand, the eigenfunction <f> a satisfies 

A(0a) = [Ci((f> a ) - C a (4> a )} - /j, a 4> a . 

Recall that from Proposition 15. 1 11 one has uniform bounds in W^' (m -1 ) on <j) a in terms 
of its L\ norm which has been fixed to 1, so that for any a G [03, 1], ||0a[|ypM( m -n < C. 
Using Proposition 13.11 and Proposition 15.111 we get 

IIA'Ml.ilm- 1 ) = C '( 1 ~ 

Using that C\ is invertible from n^L^m -1 ) to L 1 (m _1 ) we deduce that 

(5.7) ||n^ a || £ i (m -i)=0(l-a). 

We conclude the proof of (|5.5p holds for the L\ norm gathering (|5.6p and (j5.TH : 

Va6[a 3 ,l] ||0 a - 0i 1 1 2,1 < C (1 - a). 
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Let now consider the eigenfunctions & a associated to \i a for a S [«3,1] such that 
iri^a > with the normalization condition ||$ a ||wjt,i( TO -i) = 1- Proceeding similarly as 
before (by working in the space W k ' l (m~ 1 )), we can get 

[|$ a - ^iH^.i^-i) = 0(1 - a). 

Because the eigenspace associated to p a is of dimension 1, we have & a = c a (p a for some 
constant c a £ (0, oo). Then 



\ci-c a \ = \\c a a -ci<fi a \\ L i<\\$ a -$i\\ L i + \ci\\\(f)i-(f) a \\ L i=O(l-a). 

We then easily conclude that (|5.5p holds for any norm. □ 

We now use the linearized energy dissipation equation to get a second order expansion 
of the eigenvalue. 

Lemma 5.19 For a E [a3,l] ; the eigenvalues fi a satisfies (with explicit bound) 

fj, a = - p (1 - a ) + Oil - a) 2 . 

Proof of Lemma \5.19l By integrating the eigenvalue equation 

against \v\ 2 and dividing it by (1 — a), we get 

S^a) = 2 p£((f> a )- 2(1 + a)D(G a , </> a ). 



1 — a 



Using the rate of convergence of G a — > G\ and 4> a — > <f>\ established in Lemma 14.41 and 
Lemma 15.171 we deduce that 

(5.8) -^£(<Pi) = 2p£(<f )1 )-AD(G 1 ,<p 1 ) + 0(l-a). 
1 — a 

Then we compute thanks to (lA.lj) and (|A.2j) 

(5.9) £((/>!) = 2N c pel 

where cq is still the normalizing constant in (|1.33j) such that ||^i||x,i = 1. Similarly, using 
(|A.3p . (|A.4|) and the relation (jl.27p which make a link between b\ and 9\, we find 

(5.10) D(G 1 , ( / )1 ) = ^Nc p 2 el 

We conclude gathering (f5T8|) . ([5T9|) and (|5.10p . □ 
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5.6 The map a i— > G a is C 1 

The fact that the path of self-similar profiles a i— ► Gq, is C° on [03, 1] and C 1 at a = 1 was 
already proved in Lemma 14.41 Therefore we have to prove that it is C 1 for a £ [03, 1). 
Let us define the functional 

(a,g) 1 ^ ^(a,g) := Q a (g,g) - T a V v (vg). 

The map ^ is C 1 from R x (W^'^m -1 ) n C P) o) into L^m -1 ) and it is such that for 
any a 6 [ai,l), the equation ^(a,g) = has only one solution which is the profile G a . 
Moreover, for any a £ [03,1), the linearized operator D2^(a, G a ) = C a is invertible 
from W ' (m ) into L (m ) because of the spectral properties of C a established in 
Theorem 15.21 (i) & (ii) (note that here there is no eigenvalue approaching at a). Then 
using the same strategy as in Subsection 14.21 based on the implicit function theorem we 
easily conclude that a 1— ► G a is C 1 from [0:3,1) into L (mT 1 ). That ends the proof of 
Theorem ll.il (ii). 

5.7 Decay estimate on the semigroup 

We start with a lemma on non sectorial semigroups in Banach spaces. This result is a 
tool for deriving constructive decay rate on non sectorial semigroups, from the knowkedge 
on the resolvent of their generator. We do not try to prove such a decay rate for the 
semigroup in the norm of the Banach space but instead in a weaker norm (corresponding 
to the norm of the graph of some power of its generator), which shall be sufficient for our 
study of the linearized stability of the non- linear equation (|1.29[) . 

Lemma 5.20 Let A be a closed unbounded operator on a Banach space E with dense 
domain dom(A). We denote by S{t) the associated semigroup, by R(A) the associated 
resolvent set and by 1Z = 7Z(£) the resolvent operator defined on R(A). Assume that we 
have a sequence of Banach spaces E2 C E\ C Eq = E decreasing for inclusion (in most 
cases this sequence shall be provided by E^ = dom(A k ) endowed with the norm of the graph 
of A k ). We assume on the operator that: 

(i) the resolvent set R{A) contains the half plan A a for some a£l, together with the 
estimates 

sup||^(a + is)|| jBo _ > Eo < Ci 

and 

V«£l, \\K(a + is)\\ El ^E , \\n{a + is)\\E^ El < z— yr 

J- \ I S| 

for some constants C\,Ci > 0; 

(ii) the semigroup S(t) satisfies 

(5.11) Vt>0, \\S(t)\\E 2 ^Eo <C 3 e bt 

for some constants C3, b > 0. 

Then for any a' > a, there exists a constant C4 depending only on a, b, a', C\, C2, C3 
such that 

(5.12) Vt>0, \\S(t)\\ E2 ^ Eo <C 4 e a,t . 
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Proof of Lemma \5.20l We split the proof into two parts. 

Step 1. The first bound on the resolvent implies that for any x G Eq 

\\lZ(a + zs):r||£ — > 0, |s| — ► oo. 

Indeed we first consider x G dom(A) and then argue by density (since the domain dom(A) 
is dense). When x G dom(^4) the result is proved by the relation 



lZ(z)x = z 



-l 



Id + K(z) A] x. 



Step 2. Then consider the following integral of 1Z(z)x on a vertical segment with real 
part a (for some M > 0) 



ra+iM 

I M (x) := / e zt lZ(z)xdz. 
Ja—iM 



The function z — >• 7Z(z) is differentiable on this segment and we can perform an integration 
by part: 



e (a+iM)t e (a-iM)t 

I M (x) = n(a + iM)x U(a - iM)x 



a+iM zt 



— TZ(z) xdz 

a-iM * 



t t 

where we have used R'(z) = R(z) 2 . Now we estimate the Eq norm of this quantity: 



\\Im(x)\\e < 



,{a+iM) t 



t 



U{a + iM)x 



+ 



E 



,(a+iM) t 



t 



U{a + iM)x 



E 



t 



+oo 



:i + \s\ 



ds \\x\ 



E 2 - 



Therefore the integral is semi-convergent and we can pass to the limit M — > +cxd and use 
(see [32111]) that 



S(t)x 



lim 

2m M^oo 



a+iM | 

e zt TZ(z)xdz = lim Im 

_ iM 2^7^ M-»oo 



to obtain (the two boundary terms go to as M — > +oo from the first step) 



(5.13) \\S(t)x\\ EQ <C 2 — \\X\\ E2 , unu v. 



with Co = Cl 



1 



(i + H 



ds 



Using (I5TTT|) for t < 1 and (f5TT3j) for t > 1, we conclude that (I5TT21) holds with C 4 



□ 



max(C^, C 3 e b ~ a ). 

Proof of point (Hi) in Theorem \5.2[ The point (ii) of Theorem 15.21 was proved in 
Lemma ET3] and it shows that the operator C a = (Id — Uq.) C a together with the sequence 



of Banach spaces E{ 



(m 1 ), i = 0, 1, 2, for any fixed k G N and any exponential 



weight function m (as defined in (jl.28p ). satisfies the assumption (i) of Lemma 15.201 for 
any a G (fi2 , 0) . Moreover it is trivial to prove that it satisfies the assumption (ii) of 
Lemma 15.201 for some explicit b > from the decomposition C a = As — B at s(£) already 
introduced. □ 
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6 Convergence to the self-similar profile 



In this section, we consider the nonlinear rescaled equation (|1.29p and we prove the con- 
vergence of its solutions to the self-similar profile. As a preliminary let us recall some 
result on propagation and appearance of moments and regularity which is picked up from 
[24, Proposition 3.1, Theorem 3.5, Theorem 3.6]. 

Lemma 6.1 Let us consider g in G L\ n C P: q and the associated solution g G C([0, oo); L3) 
to the rescaled equation U.29\) . Then 

(i) For any exponential moment weight m (as defined in \1.28\) ) with exponent s G 
(0,1/2) and any time to £ (0>°o); there exists a constant Mi = Mi (to) such that 

(6.1) sup HffC^OIU^m- 1 ) < M i- 

[to, 00) 

Moreover, if g m G L 1 (m _1 ) /or some polynomial or exponential (with exponent s G 
(0, 1) ) moment weight m then \6. 1}) holds (for this weight m ) with to = and some 
constant M\ = M 1 (\\g in \\ L i( m -i)). 

For the two following points we now assume that for some constants c\,T G (0, 00) 
there holds 

(6.2) inf > C1 , 

and we state some smoothness properties of the solution g which depend on c\ but 
not on T nor a. 

(ii) Assume $6.2\) . Then for any ko G N there is qo = qo(ko) G N such that «/||s , m||_f/ fc onL 1 — 
Co holds, then for any c\ G (0, 00) there exists C\ = C\(Cq,ci) G (0,oo) such that 
for any time T G (0,oo), we have 

(6.3) VtG[0,T], ||<?(V)llif*i <Cl, 
with ki = if k = and k\ = k — 1 if k G N* . 

(Hi) Assume and that g m G L 2 , with WgmW^nL 1 — Mi G (0, 00). Then there exists 

A G (—oo,0) and for any exponential weight function m with exponent s G (0, 1/2) 
and any k G N, there exists a constant K (which depends on p,c\,M\,k,m) such 
that we may split g = g s + g R with 

(6.4) Vt G [0,T], H/CMIUolH™- 1 ) ^ K > \\9 R (tr)hl < Ke xt . 

Remark 6.2 It is worth mentioning that these estimates are uniform with respect to the 
inelasticity parameter a G (0,1). Indeed, one the one hand, this was already the case for 
the moment estimate h6.1\) in \2lj\ Proposition 3.1]. On the other hand \6. 3\) and {6.$ 
from f2~4\ Theorem 3.5, Theorem 3.6] were (partially) based on the use of the damping 
effect of the anti-drift term (whose coefficient was fixed to r = 1). Here the damping effect 
of the anti-drift term vanishes (r a — > 0) but it is replaced (as for the elastic Boltzmann 
equation) by the lower bound on the energy h6.2\) which allows for a control from below 
on the convolution term L(g) appearing in the loss term of the collision operator (see 
Lemma \2.3\) . which is enough to conclude also in this case. 
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6.1 Local linearized asymptotic stability 

Let us first consider the nonlinear evolution equation (|1.15p in Pi H k , and the 

associated equation on the fluctuation h of a solution g around the unique equilibrium 
G a : g = G a + h and 

d t h = C a h + Q a (h, h). 
Let us start by stating an inequality that we shall need in the sequel. 



Lemma 6.3 For any exponential weight function m (as defined in \1.28\) ), there is a 
constant C G (0, oo) such that for any h G W^' l (m~ l ) and any a G (0, 1), 

||n a Q a (/i,/i)|| L i (m -i) < C(l -a) IWI^i^-i)- 

Proof of Lemma 16.31 We write 

n a Q Q (h, h) = U a (Q a (h, h) - Qi(/t, h)) + (n Q - ILjQ^/i, h). 
On the one hand, from Lemma 15.151 (i) and (|3.4p . there is C G (0, oo) such that 
\\n a (Q a (h, h) - Q X (H, h))\\ L i {m -x } < C (1 - a)\\h\\ 2 w s, 1{m _ iy 

On the other hand, from ()5.3[) and (|3.ip . we get 



||(n Q - nOQiCfc, /»)|Ui (m -i) < c(i - a) \\h\\ 2 w s, 1{m _ iy 

The proof of the lemma is immediate by gathering the two previous estimates. □ 
We now state a first local linearized stability result. 

Proposition 6.4 For any a G [a^, 1), the self-similar profile G a is locally asymptotically 
stable, with domain of stability uniform according to a £ [a^, 1). 

More precisely, let us fix p G (0, oo) and some exponential weight function m as 
in il.28\) . There is k%,q\ G N* such that for any Mq G (0, oo) there exists C, e G (0,oo) 
such that for any a G [03, 1], for any g m G H kl n L l (m~ qi ) with mass p, momentum 
satisfying 

(6-5) Il5in||j2-*in£ 1 (m-«i) — «9in ~ G a ||_Li( m -i) < £, 

the solution g to the rescaled equation U.29\) with initial datum g m satisfies 



(6-6) W>0, ^(gt-G^h^^KCy^-Gahi^e^, 

L (m- 1 ) — C I bin — GaWv-imr 1 ) 



(6.7) Vt > 0, ||(Id - n a ) (g t - G^IUm- 1 ) < C \\9in ~ Gahi^-i) e^^ 1 . 



Proof of Proposition ROl Step 1. Let us first denote by c\ the constant given in Step 5 
of Proposition 12.11 such that 

V«e[ai,l), £(G a )>2 Cl . 
We may then fix £q G (0, 00) in such a way that 
(6.8) \\g - G a || L i( m -i) < e implies £(g) > a, 
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and define 

% := sup{T, £(g t ) > c x Vt G [0,7]} G (0,oo]. 

From Lemma IBTTI (i) & (ii), there exists M G (0, oo) (depending on p,ci,ki,q±, Mq) such 
that for any T G (0, oo) there holds 

(6-9) sup \\g\\ H k inLHm - qi) < M. 

te[o,T»] 

Let us now consider the fluctuation ht — Qt~ G&. Thanks to the mass and momentum 
conservations, it satisfies h t G Co,o for all times, as well as the bound (|6.9p . We define the 
following decomposition on h: 

h 1 = U a h and h 2 = (Id - IL a )h =: IL^h. 

Since the spectral projection n a commutes with the linearized operator C a , the equation 
on h 1 writes 

d t h l = p a h l + U a Q a (h, h). 

Multiplying that equation by (sign/i)m _1 and integrating in the velocity variable, we 
deduce thanks to Lemma RT31 and to (1B.2|) . (|6.9|) that on (0,T*) the following holds 

^ll^lli^m- 1 ) < Ma II h l \\ui(m-i) + C (1 -a) 11^11^-3,1 (m _i) 

(6.10) < (l-ajJCill^ll^ + C-ill^H^-Call^lUi^-i) 

for some constants C\ depending on M and the possible choice C2 = p/2 for Ci- For the 
second part h 2 we have the following equation 

d t h 2 =Iij; C a h 2 + U^Q a (h,h). 

Since the linearized operator C a restricted to 11^ generates the semigroup R a (t) defined 
in point (iii) of Theorem 15.21 the Duhamel formula reads 

h 2 (t) = R a (t)h in + [ R a (t- s)U^Q a (h,h)(s)ds. 
Jo 

From (15. ip and ( 13.11 ) we have 
\\h 2 



LHm-i) <Ce^\\h- m \\ LHm -, )+ C f\^ t ^\\h(s)\\ 2 wr[m _ 1) ds. 

J 



We deduce 



(6.11) ||/ i ?||Li( m -i)<C 3 e^||^nll+C4 J ^^(KftV^ + ll^lliV 1 ))^ 

with C4 depending on M thanks to (|B.2[) and (|6.9p . It is then easy to show by comparison 
arguments from (|6.10p and (|6,lip that there are < £2 < £\ < £q (one can take for 
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1 /2 1 /2 

instance e\ < Eq/2 satisfying 2 C% e x < C2 and 2 C4 £-/ < 1/2 and next £2 < £1 satisfying 
C3 £2 < ^i/2) such that 

H^hJIzHm- 1 ) + ll^fnllL^m- 1 ) < ^2 implies sup maxlll^Hii^-i), ||^||ii( m -i) } < £i- 

te[o,T»] L J 

(6.12) 

Gathering (|6.8p and (|6.12p we deduce that there exists e E (0, £2) such that under condition 
(I6.5P there holds T* = 00 as well as 

sup \\g- Gy z i (m -i) < 2ei < e . 
te(o,oo) 

Step 2. In a second step, coming back to (|6.1ip and to the integral version of (|6.10p and 
setting y(t) = \\h \\ + \fi a \ \\h 2 \\, we obtain 

(6.13) y(t) < C 5 1 y(0) + C 6 | Ma | / e*« ^ y{sf' 2 ds. 

Jo 

Then we have to the following variant of the Gronwall lemma whose proof is the same 
that the one of [271 Lemma 4.5] and is therefore skipped: 

Lemma 6.5 Let y = y(t) be a nonnegative continuous function on M + such that for some 
constants a, b, 9, fj, > 0, 



y(t) < ae-^X + b (J e^ ( *- s) y(s) 1+e ds) 



(as compared to \2T\ Lemma 4-5], X needs not necessarily be y(0)). Then if X and b are 
small enough, we have 

y(t) < CXe~^. 

for some explicit constant C > 0. 

Thanks to the uniform smallness estimate on y(t) we can apply the lemma with 9 = 1/4: 
for instance, and we get 

y(t)<CVy(0)e^* 

from which we deduce the estimate (j6.6[) for the h 1 part of g — G a . Finally, we may insert 
that estimate on h 1 in (|6.1ip and we get 

ll>» 2 (f)[Ui(m-i) < CUe^ + e^^) ||^|Ui (m -i) + C, \\h\s)fl 2 {m _ 1) ds. 

The same kind of computation yields to 

\\h 2 \\ < CaeW)"** \\h(0)\\ 
from which (16.71) follows. □ 
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6.2 Nonlinear stability estimates 

In this subsection we shall prove that when the inelasticity is small, depending on the size of 
the initial datum (but not necessary close to the self-similar profile), the equation (|1.15[) 
is stable. This relies mainly on the fact that the entropy production timescale is of a 
different order (much faster) that the energy dissipation timescale as a — * 1. This point 
is familiar to physicists (see for instance [9]) which separate, for granular gases with small 
inelasticity, the molecular timescale (the level where entropy production effects dominate) 
and the cooling timescale (much slower than the molecular timescale). 

Proposition 6.6 Define £?2 := max{A;o, k\}, qi := maxjqo, Qi, 3}, where ki and qi are 
defined in Theorem 13.51 and Corollary \3.4\ For any p, £q, Mq there exists 04 G [03,1), 
ci G (0,oo) and for any a G [a$, 1] there exist tp = p(a) with p(a) — > as a — > 1 
and T = T(ot) (possibly blowing-up as a — > 1 ) such that any initial datum < g in G 

L\ 2 n H k * n C pA s with 

llffinlLi 2 ni? fc 2 < M), 

i/ie solution g associated to the rescaled equation tl.29\) satisfies 

Vt>0, £{g t )>c 1 
and for all a' G [04, 1) and then all a G [a 1 , 1] 

(6.14) Vt>T(a'), Hst-GaLi^cO. 

Proof of Proposition 16.61 Let us consider a solution 5 G C([0, 00); L^ 2 n // fc2 ) to the 
rescaled equation (|1.29j) with given initial datum g m , whose existence has been established 
in [231 121] • We split the proof of the Proposition into five steps. 

Step 1. From the propagation and appearance of uniform moment bounds [231 Proposi- 
tion 3.1, (hi)], which it is worth noticing have been obtained uniformly with respect to the 
elastic coefficient (see also [8]), there exists C\ G (0, 00) such that 

(6.15) sup||g|| L i < Ci. 
Let us define c\ := min{£ (Gi), Sq}/4, and 

(6.16) % :=sup{T ; Vi G [0,T], £(g(t,-)) > c x }. 
Next from the equation on the evolution of energy 

(6.17) S'(t) = -(1 - a 2 ) 61 D e {g) + {l-a)2pS 
and ()6.15p there holds 

\S'(t)\ < C 2 (l-a) Vt>0 
(take for instance C2 = 2b\ C 2 + Ci), from which we deduce that we necessarily have 

T* > C 3 (1 - a)" 1 

(take for instance C3 = (3/4) E§jC<i). 
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Step 2. From point (ii) of Lemma 16. 11 we have for some constant C5 G (0, 00) 

(6.18) ViG[0,T*] \\g t \\ Hk2 <C 4 . 

Moreover from Lemma 12.6} for any time t\ £ (0,T*), there exists some constant C5 = 
C§{p,Ci,t{) such that 

(6.19) Vie fa, T*] ff (t,v) > C^e^H 8 . 

S'tep 5. With the notations of Theorem 13.51 we compute the evolution of the relative 
entropy of g(t, •) with respect to the associated Maxwellian M[g(t, •)], and we obtain 



±H i9 \M [9]) = J t m-±l N ^M [9] =*H {9) 



2£ dt 



p N 

= -D H:a (g)-tj(l-a)D e ( 9 ). 
Next from Lemma 13.41 and the estimates (|6.15p , (|6.16p , (|6.18p and (|6.19p we have 

j t H(g\M[g\) = -D Hil (g,g) + 0(1 - a) on (t 1} T.). 

Then from (|3.1ip . we are then led to the following differential inequation on the relative 
entropy 

j t H(g\M[g\) < -C 6 H(g\M[g}) 2 + C 7 (1 - a) on (t u T*). 

By straightforward computations we deduce that independently of the value ofH^g^ 

(this "loss of memory" effect is typical of differential equations with overlinear damping 

terms), we have 



Vt€[ti,T„], H(g t \M[g t ])<C 8 (l-a) 



1/2 l + e -C 9 (l-a)V 2 (*-M 

l _ e -c 9 (i- a y/ 2 (t-H) 



for some explicit constants. As a conclusion, defining £2 := t\ + C Q 1 (1 — a) l l 2 and 
choosing a G [03, 1) in such a way that t-i < T* we have for a G [a', 1) 



VtGfa,T*] ^( 5 (t)|M[g])<Ci (l-a) 



1/2 



Finally, using Csiszar-Kullback-Pinsker inequality (|3.10p . as well as Holder inequality, we 
obtain under the same conditions on a and the time variable: 

(6.20) ||s - M[g)\\ Ll <C\\g- M[g]\\ l ff \\gfg < CH(g\M[g]) 1 ^ < 0(1 - a) 1 / 8 . 

Step 4- Now let us go back to the energy equation (|6.17p . First, with the help of the 
moment bound (|6.15p . one may write 

£'(t) = 2(1 -a)[p£-b 1 D £ {g) + 0(1 - a)]. 

Thanks to (|6.20p we deduce 

£'(t) = 2{l-a){p£-b 1 D e (M[g]) + 0{(l-a) l l*)) on (i 2 ,T*). 
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Finally, thanks to (|3.16p . (|3.17p and the relation £{g) = pN6(g), we get on (t2,T*) 

(6.21) S'{t) := E(£(t), a) = (1 - a) [k 3 £ {£{ /2 - £ 1 ' 2 ) + 0((1 - a) 1 / 8 )], 

where £\ = pN Q\ with #1 is the quasi-elastic self-similar temperature defined in (|1.27p . 
We may then choose a" G [a', 1) such that E(ci, a) > for any a G [a", 1). We conclude by 
maximum principle that T* = oo for a G [a", 1). In particular, all the previous estimates 
on g are uniform on fa, oo). 
•Step 5. Thanks to (|6.21j) we easily get 

j t {£- £i? < -(1 - a) [k 5 {£ - £,? + - af/% 
so that (for some constants a, b > 0) 

Vt>t 2 , \£{t)-£i\ < \£{t 2 ) - £i\e- a ^- a){t - t2) + b{\ - af l% . 
Setting T(a) = max{t2, c (1 — a)" 1 } for some suitable constant c > 0, we then obtain 

(6.22) \£~£i\ = 0{{l -a) 1/8 ) on [T(a), oo). 
In order to conclude that (|6,14p holds, we write 

g(t) — G a = (g(t) - M[g(t)}) + (M[g(t)] - G x ) + {G x - G a ), 

and we estimate the first term thanks to (|6.20p . the second term thanks to (|6.22p and the 
third term by ^TFS). □ 

6.3 Decomposition and Liapunov functional for smooth initial datum 

The proof of the gobal convergence (point (v) of Theorem II. lh for smooth initial data only 
amounts to connect the two previous results of Propositions l6.4l and l6.6l bv choosing a such 
that ip(a) < e where e is the size of the attraction domain in Proposition 16.41 and <p(a) 
is defined in Propositions 16.61 More precisely, we state without proof the straightforward 
combination of Propositions 16.61 and Proposition 16.41 

Corollary 6.7 Let us fix an exponential weight function m as in il.28\) . with exponent s G 
(0, 1). Then for any p, £q, Mq there exists C and 05 G [0:4, 1) (depending on p, £q, Mq,iti) 
such that for any a G [«5, 1) and any initial datum < g m G L l (m~ q2 ) n H k2 satisfying 

9 m G C p 

,0,£o' 1 1 Sin Hi 1 {m~ q 2 )nH k 2 

<M , 

the solution g associated to the rescaled equation U.29\) satisfies 

Vt>0, IIIkGfc-GeOH^-i) <Ce^\ 

Vt>0, || (Id - H«) (g t - G a )|| Ll(m -i } < Ce^^K 

Remark 6.8 Note that the constant C in the rate of decay does not depend on a. This 
comes from the fact the size of the linearized stability domain is uniform as a goes to 1 
in Proposition \6.4\ which allows in Propositions \ 6.6\ to pick a fixed a' such that in the 
estimate {6. 14\ ) tp(&') is less than this size, and therefore that the time T(a') required to 
enter this neighborhood does not blow-up as a goes to 1. 



63 



As a by-product of the previous propositions, we state and prove a result which pro- 
vides a partial answer to the question (important from the physical viewpoint) of finding 
Liapunov functionals for this particles system. Let us define the required objects. We 
consider a fixed mass p and some restitution coefficient a whose range will be specified 
below. At initial times, non-linear effects dominate and therefore we define 

Hi (g) :=H(g\M[g]) + (£-£ a ) 2 

where £ a = £{G a ) is the energy of the self-similar profile corresponding to a and the mass 
p. At eventual times, linearized effects dominate. Therefore we define a quite natural 
candidate from the spectral study: 

n 2 (g) ■■= H^Hii^-i) + (!-«)/ \\Ra(s) h 2 \\ L2 ds, 

Jo 

with h 1 = H a h, h 2 = h and h = g — G a . 

Proposition 6.9 There is k^ G N big enough (this value is specified in the proof) such 
that for any exponential weight function m as defined in lll.28\) . any time to G (0, oo) 
and any p,£q,Mq G (0,oo), there exists G (0, oo) and a§ G [as,l) such that for any 
a G [a^, 1] and initial datum g m G H ki n L 1 (m _1 ) satisfying 

gin G C p ,o,£ , Ibinll^nLHm- 1 ) ^ M 0, <7in(«) > Mq 1 e~ M ° \< , 

the solution g to the rescaled equation il.29\) with initial datum g m is such that the func- 
tional 

H{g t )=n 1 {g t )l {Hi{gt) ^ ] +n 2 {g t )l {Hi{gt) ^ ] 

is decreasing for all times t G [0, +oo). Moreover, TC(g(t, •)) is strictly decreasing as long 
as g(t, ■) has not reached the self-similar state G a . 

Proof of Proposition 16.91 We split the proof into three steps. 

Step 1: Initial times. Taking k^ > k 2 and a G [04, f ), we know from the proof of Proposi- 
tion HOI] that the solution g satisfies that 

Vt€[to,oo), Hs(t,0llH*4nLi(m-i) <Mi, g(t,v)>M^e- M ^ 8 , 

for some constant M\ G (0, 00) (recall that 04 was adjusted in terms of p, £0, Mq). Coming 
back then to Steps 3 and 4 in the proof of Proposition I6.6L we obtain the two following 
differential equation on (to, 00) 

j t H(g\M[g]) < -K x H(g\M[g}) 2 + 0(1 - a) 

and 

j£ = 2p (1 - a) [K 2 £ (£ l J 2 - £ 1 ' 2 ) (£ - £ a ) + 0((1 - a) 1 / 8 ) 

for some constants K{ G (0, 00). We easily deduce that for any k G (0, 00) there exists 
a K G [05,00) such that 

(6.23) ^t^ 1 ^ < for any 1 £ (°' 0O ' ) SUch that Hl ^ 9i} > - K - 
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Step 2: Eventual times. Let us first remark that from point (iii) in Theorem 15.21 (iii) and 
the interpolation inequality (|B,2j) . for any q G N* there exists k,k' G N and Cj G (0, oo) 
such that 

||-Ra^ ||^2 — Cl ||-Rq^ 1 1 ^fc,i (m,— <?/2) 

< C 2 e^ s \\h 2 \\ w k+2,i {m _ q/2) < C 3 e^ s \\h\\ Hk ' nL i {m - q) , 

so that, taking k^ big enough, the functional Ti.2(g(t, .)) is well-defined for any times 
t G (0, oo). First observe that from (|6.10p there holds 

(6.24) J^WlHm-l) <(!"«) [*1 ll^llKV 1 ) " ^ 2 ll^lli^m- 1 ) 

Second, we compute (with the notation of Subsection 15.71) 

-/ 11^(5)^11^^ = 2/ / (e 8 ^ /i 2 ) [e a£ -(A»fc 2 + niQ«(/i,/i))]<fa du- 
al Jo Jo Jr n 

On the one hand, 

+oo n 

/ (e sCa h 2 ) [e sCa C a h 2 ]dsdv 

On the other hand, 

"+0O 



r+oo i- 

h = 2 / / (i? a ( s )/i 2 )[i? a ( s )n^Q a (/ i ,/ i ))]d s ^ 

/■•+•£» 

< 2C, iy o ll i? a( S ) /l2 |liy fc i. 1 (m-' J '' 2 ) ll i? a( S ) n a ( 3Q( /l ' /l )llTyfei.i(m-'j/2)ds 

< C^Qf e 2 ^ s ds^ ||/i 2 ||^ 1+14(OT _ 9/2) ||Qa(^^)||^i+i.i (ni _ 9/2) 

< C3 H/i 2 !!^ 4 ll^ 2 H^t nL i (m -i) INI?? H^llfffcaniiCm- 1 )' 

for some k% G N given by Proposition IB.ll Taking k^ > ^3, we then obtain 

(6.25) I II^W^H 2 ds < K 3 \\h\\f - \\h 2 \\\ 2 . 



dt 



Gathering (|6.24p and (|6.25p and using some interpolation again, we deduce that there 
exists k! G (0, 00) such that 

(6.26) -^T~(-2(gt) < for any t G (0, 00) such that ||fot||.Li < k'. 

Step 3. We conclude putting together (|6.23p and (|6.26p . and using (|3.10p . (|3.12p in order 
to prove that 

T~L\(g) < K implies H/hIIl 1 < 
for a G [«6) 1] for some a$ G [as, 1). □ 
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6.4 Global stability for general initial datum 



We first state and prove a regularity result on the iterated gain term which is the inelastic 
collision operator version of the same result proved for the elastic collision operator in 

nan]. 



Lemma 6.10 There exists a constant C such that for any f,g,h G Ll(R N ) and any 
a G (0, 1] there holds 



(6.27) 



\Qi(f,Qi(g,h))\\ L3 <C\\f\U\\g\ 



Li' 



Proof of Lemma \6.1(K We follow [261 lemma 2.1] and [H lemma 2.1] and we make 
use of the Carleman representation introduced in |24^ Proposition 1.5]. Let us consider 
f,g,h G L^M^) and <fi G L°°(M W ). We apply twice the weak formulation of the gain term 



Q + (f,Q + (g,h))(v)<p(v)dv 



Q + (g,h)(v) 



g(v)h(vi) f(v 2 ) 



f(v 2 ) \v - v 2 \ 



s 2 



(w' 2 ) da 2 dv 2 



\v — V\\ 



3- 



\v 1 - v 2 \ 



dv 

<{v 2 ) do 2 do-\ 



dv dv\ dv 2 



with w 2 = V'(v 2 ,v,cr 2 ), v[ = V'(v,vi,ai) and therefore v 2 = V'(v 2 ,v[,a 2 ). Recall that 
for any given v,v*,a G Mr, we define 



w 



V + V* 



1 + e 



u = v — V* 



(1 — 7) u + 7 \u\ a 



and then 



7 



V' = V'(v, v m , a) = - + — = v + ± a-u) 

V* = V*{v,v*,a) = — - — = v*-- {\u\a-u). 

We denote by $ = $>(v, vi,v 2 ) the term between brackets in the last integral. Introducing 
the point w\ and the set S VjVl>E defined by 



w\ := (1- 7/2) v + (7/2)^1, Sv !V1;£ 
we get 



z G 



z — w\\ — (7/2) \v — Vl\ 



< e/2} , 



(6.28) $ 



lim 



\V — V\ \ £^0 £ 



R N JS 2 



1 S v , vl , e (v[) K - v 2 \ c/)(v' 2 ) da 2 dv[. 



Remarking that v 2 = v 2 + (7/2) (\u 2 \ o~ 2 — u 2 ) with u 2 = v[ — v 2 , we observe that the 
integral term \& e is very similar to the collision term Q + (here v 2 (resp. v\, a 2 , 7, v 2 ) 
plays the role of v (resp. v\, a, (3, 'v) in the gain term) and therefore we may give a 
Carleman representation of \P £ , The same computations as performed in [24|. Proposition 
1.5] yield 

V F = — I I 1Sv, vl ,e( v 'l) \ v 2 ~ ^l" 1 4>{v 2 ) dE{v'l) dv 2 



T 



v 2 ,v 2 



66 



where E V2 y/ is the hyperplan orthogonal to the vector V2 — v 2 and passing through the 
point £l V2jV » = V2 + (1 — 7 _1 ) (v 2 — v 2 )• Here v'l stands for the post collision velocity issued 
from v[, that is v'% = V*(y2,v' 1 ,o~2), and then, thanks to the momentum conservation, 
v i := v 2 + v 3 ~ v z- We finally define EL^y the hyperplan orthogonal to the vector V2 — v' 2 ' 
and passing through the point O' „ = v 2 + (1 — 7 _1 ) (^2 — v 2 ) and we get 

(6.29) * e = 4 / / l^^OI^'-^rVK)^!)^. 

Now, arguing as in [U lemma 2.1], we see that the measure of the intersection T £ between 
the plane It V2 y> and the thickened sphere S V)Vl>£ is bounded by ire~/\v — v±\ and that 
v 1 G T e implies that v 2 G B £ with 

5 s := jz G M^; |z| 2 < |t>| 2 + \vi\ 2 + 2e(\v\ + \v x \) + e 2 1 1^2 1 2 } - 

Gathering these estimates with (|6.28p and (|6.29p we get 



V — V\\ £->0 £ J r n \V 2 — V2\ 

7 d e-^Oj^N \v'{ ~V 2 \ 7 d Jjjiv |t# - «a| 

where we have defined B° := {z G K^; |z| 2 < |w| 2 + |fi| 2 }. Using [TJ lemma 2.2] we may 
conclude as in the end of [U lemma 2.1] and therefore (|6.27|) follows. □ 

We second establish that the solution g of the rescaled equation (11.290 decomposes 
between a regular part and a small remaining part as it has been proved for the elas- 
tic Boltzmann equation in [28], and then partially extended to the inelastic Boltzmann 
equation in [24] . As compared to this last paper, this result relaxes the assumption on 
the initial datum to g m G L\, but at the price of the hypothesis of a lower bound on the 
energy. 

Lemma 6.11 Consider g m G L\ and the associated solution g G C([0, 00); L\) to the 
rescaled equation 111.29]) . Assume that for some constant p,c±, M±,T G (0, 00) there holds 

(6.30) g in G C Pi0 , WginWq < Mi, Vt G [0,T], £(g(t,-)) > Cl . 

Then, there are 017 G [«6>1) an d A G (— 00, 0), and for any exponential weight function 
m (as defined in il.28]) and any k G N, there exists a constant K (which depends on 
p, ci, Mi, k, m) such that for any a G [07, 1], we may split g = g s + g R with 

(6.31) VtG [0,T], ||/(t,-)||ff*nLi(m-i) < K, \\g R (t, .)\\ L i < Ke xt . 

Proof of Lemma 16.111 The starting point is to write the rescaled equation (|1.29p in the 
following way 

dg 

— + T a v V v g + £g = Q+(g,g), 
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with £(t, v) := T a N + L(g(t, -))(v). Introducing the linear semigroup 

(U t h) (v) = h{v e' Ta * ) exp (- J £(s, v) ds 

and using the Duhamel formula, we have 

9t = U t g- m + \ U t - s Qt(9s,g s )ds. 
Jo 

We iterate that last identity and we obtain g = g± + gf with 

r-t pt rs 

9l=U t g- m + Ut-sQi(g s ,U s g in )ds, gf = U t - s Q+(g s ,U s - u Q+(g u ,g u ))duds. 
Jo Jo Jo 

On the one hand, the energy lower bound (|6.30p and Lemma 12.31 imply that there exists 
a constant C2 G (0, oo) such that 

{U t h)(v) < e~ C2t (V it h)(v) with (V^h)(v) = h[£v) and & = e~ Tat 

On the other hand, straightforward homogeneity arguments leads to 

Qt(g, Vsh) = r N -' Vt-i Qi(v^ g , h) 

and |.| 9 ||lp = £~ q ~ N /P ||/t|.| 9 ||ip for any functions g, h and positive real £. From these 
considerations we deduce that 

\\g?(t)\\v < e^'^ \\g ia \\ L1 + e ((^+D— >)« || 5in || XJ sup || 5s || x ; < C e~^\ 

1 s>0 

for some constant C and for any (1 — a) small enough. In the same way, we have 

lbfWllL3 < /"*r e [(2iV/3+1} ^- C2l( *- CT) r x g s ,Qi(g„,g a ))\\ L 3dads. 

Jo Jo S " CT 

Taking (1 — a) smaller if necessary and using Lemma 16.101 we obtain 
lbf(i)llL3 < ( f f e-^ /2){t - a) dcjds\^ V \\g s \\l,, 

\JoJo J s>0 2 

which ends the proof of (|6.3ip in the case k = 0, with the help of point (i) of Lemma 16. 11 
The general case k £ N* is then treated by following the strategy introduced in [28J and 
using the result of appearance of regularity proved in [24] (and recalled in point (iii) of 
Lemma 16. ip . □ 

We third recall a classical L 1 stability result for the elastic Boltzmann equation which 
has been established in [24, Proposition 3.2] for the rescaled equation (|l,29p . 

Lemma 6.12 Consider < gl n ,gf n 6 i^fl C p ,o and the two associated solution gi 6 
C([0, oo);L\) n L°°(0, oo; L\) to the rescaled equation hl.29\) . There exists C sta b € (0, oo) 
(only depending on b and sup 4>0 \\g x + g 2 ]^^ such that 

Vt > 0, || 9 2 - gl\\ L x < || ft 2 n - gi\\ L i e c stab*. 
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Proof of point (iv) of Theorem \l.l[ Let us consider g m E L\ n Cp,o,£ in with H^mlLi < Mo 
for some fixed <?i n , Mo E (0, oo) and g the associated solution to the rescaled equation 
(|1.29p which has been built in [23]. We know that there exists Mi E (0, oo) such that 

(6.32) sup \\g{t,-)\\ Ll <Mi. 

(0,oo) 3 



Step 1. We define 

T* := sup{T E (0, oo), £(g(t, •)) > ci Vt E [0,T]}, ci := min{£ in , £ i}/2. 

We shall prove that T* = +oo. We argue by contradiction, assuming that T* < oo. From 
the equation on the energy (|6.17p and the uniform estimate (|6.32p and from the definition 
of T* we have 

(6.33) % > d (1 - a)' 1 and £'{%) < 0. 
Thanks to Lemma 16. Ill we may decompose 

g = g s + g R on (o,ti), 

with ti E (0, T*) to be fixed. At time ti we initiate a new flow starting from the smooth 
part of g. More precisely, we decompose 

g = g s + g R on (ti,r„), 

with g s (t\) = [p/p(g (tx))]g (ti), g s solution (with mass p!) to the equation (|1.29|) on 
(ti,T*) and g R := g — g s . On the one hand, from (|6,3ip and Lemma 16.121 we have 

\\9 R (t)\\ L i < Ce c stab(^-*i)+^i on (ti,T*). 

We choose t\ = rjT^ with E (0, 1) in such a way that C sta b (1 — rj) + X r\ = A/2. We have 
then proved 

(6.34) ll^(i)]ILi <^e (A/2)Cl(1 " Qrl on (i x ,T*). 

On the other hand, following Step 3 in the proof of Proposition 16.61 we deduce a similar 
estimate as (|6.20p . namely 



(6.35) ||5 S m, •)- M[g s (T«, -)]|| L i = C?((l - a 



for any (1 — q) small enough chosen in such a way that the intermediate time £2 defined 
in Step 3 of the proof of Proposition 16.61 satisfies t\ + £2 < 7*. Gathering (|6.34p and (I6.35|) 
we obtain 

\\~g(T„-)-M[g(n,.)}\\ q =0((l-a)^). 

Coming back to the equation (|6.17p on the energy and proceeding like in Step 4 in the 
proof of Proposition 16.61 we get 



£'(%) > (1 - q) h Cl (£l /2 - c{ /2 ) - C (1 - a) 



> 



for any (1 — a) small enough. That is in contradiction with (I6.33P and we conclude that 
% = +00. 
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Step 2. Thanks to the previous step, we have a uniform in time lower bound on the 
energy, and therefore we can run the decomposition theorem for all times. 

By applying the decomposition theorem as in Step 1 for a given time t £ (0, do), 
starting a new flow at t\ = rjt taking [p/ p(g s (t±))] g s (t\, •) as initial datum, and then 
using Corollary 16.71 on the smooth part g s (s,-), s G [ti,t], we find that at time t, the 
solution gt decomposes as gf + g^, where approaches the self-similar profile with rate 
q e p. a (t-ti) ) foofc [ s q e (i-r]) fiat ^ anc j ^he remaining part g R goes to with rate Ce^ 2 '*. 
Since |A/2| is larger than (1 — rj) \p a \ for (1 — &) small enough, it concludes the proof of 

CGEj) . □ 
A Appendix: Moments of Gaussians 

We state here some moments of tensor product of Gaussians. 
Lemma A.l The following identities hold 

(A.l) f Mi o i \v\ 2 dv = N, 
Jr n 

(A.2) f Mi oi \v\ A dv = N (N + 2), 

JR N 

(A.3) f M 1A1 (M 1A1 )*\u\ 3 dvdv* = 2 3 / 2 [ M 1A1 \v\ 3 dv, 
Jr n xr n Jr n 

(A.4) / Mi, ,i (Mi, ,i)* M 2 \u\ 3 dvdv* = V2(2N + 3) [ M 1A1 (v) \v\ 3 dv . 

JR N xR N JR N 

Proof of Lemma I A. 1 I The proof of (|A.1|) and (|A.2|) being straightforward and the proof 
of (|A.3p being very similar to the proof of (|A.4|) we only prove (|A.4|) . We first notice that 



Ml,o,i (Ml,o,i)* M 2 \u\ 3 dv dv* 

xR N 

1 r 

Mi, 0l i (Mi )0 ,i)* (M 2 + M 2 ) \ufdvdv* 

lR N xR N 

Mi )0 ,i (Mi )0 ,i)* (\v + v*\ 2 + \v — v*\ 2 ) \u\ 3 dv dv*. 



R^xR^ 



Making use of the change of variable (v,v*) — > (x = (v + v^)/y/2, y = (v — v # )/y/2), we 
then get 



/ Mi (Q) i (Mi i0 ,i)* |f| 2 |u| 3 cfocfc* 

JR^xR^ 



V2 / M lj0il (x)M h0tl (y)(\x\ 2 + \y\ 2 )\y\ 3 dxdy 

JR N xR N 

V2N / Mi. ,i(v)|w| 3 dt; + v / 2 / Afi j0 ,i(t?) M 5 cfo 
Jr n Jr n 

V2(2N + 3) / Mioi(w) \v\ 3 dv. 
Jrn 

□ 
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B Appendix: Interpolation inequalities 



Lemma B.l (i) For any k,k*,q,q* G Z with k > k* , q > q* and any 9 G (0, 1) there 
is C G (0, oo) such that for h G W q »* ' (m _1 ) 

(B.l) \\h\L,k,i, ix <C||/i|| l - l ?. 1 \\h\\ e k ,* lr „. 

with k**,q** G Z sucfc ffc<rf fc = (1 - 9) k* + 9 k** , q = (1 - 0) g* + 0g**. 

(ii) For any k,q G iV* and any exponential weight function m as defined in \1.28\) . there 
exists C G (0, oo) such that for any h G FL k% n L^m -12 ) with fc* := 8A: + 7(1 + AT/2) 



( B - 2 ) IWIw^fm-^ - C 



1/4 ii i ii 1/4 ||7||3/4 



Wf'^m- 1 ) ~ w H'"!!^ H^llLitm-i 2 ) ll'^llLiCm- 1 )- 

Proof of Lemma lB.ll The inequality (|B.1|) in point (i) is a classical result from interpolation 
theory. Let us focus on point (ii). We prove the inequality (|B.2[) for h G <S(R ) and then 
argue by density. On the one hand, we observe that for any t there exists C such that 

\\h\\ 2 Hl < C \\h\\ L i \\h\\ H # , tf:=2l+l + N/2. 

Iterating twice that inequality, we get (for some related exponents k\k^) 



(B.3) \\h\\ Hk , <C\\h\\f\\h\\f kt . 



On the other hand, using first Cauchy-Schwartz inequality, plus the same argument as 
above and Holder's inequality, we obtain 

|V2 ||, ||l/2 



(B.4) < C\ lhu1 ^ 8 \\h\\ 3 / 8 ii aii 1 / 2 



We conclude gathering (jB.4[) and (jB.3|) . □ 
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